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ABSTRACT 


This report describes the initial phases of an investigation of hybrid 
propellant combustion instability. Experimental studies carried out during 
the first year of this program were principally concerned with the delinea- 
tion of steady-state hybrid propellant regression rate /pressure coupling 
in the pressure-sensitive regime, because acoustic instability is expected 
to depend on the same coupling mechanism. The data established that the 
onset of the observed dependence of regression rate on pressure at low 
pressure is attributable^to the behavior of the chemical kinetic processes 
in the gas-phase flame zone. This result provided the basis for theoreti- 
cal studies performed during the program, 

A theoretical model, based on classical turbulent flame theory, was 
proposed and a mathematical analysis was developed that provided good agree' 
ment with the observed steady-state regression rate/pressure dependence at 
the fixed oxidizer mass flux. 

During a few of the experiments, a spontaneous instability correspond- 
ing to the longitudinal mode of the chamber was observed to develop. Fu- 
ture studies will be concerned with the delineation of conditions under 
which the regression rate response to pressure fluctuations can support 
such instability. 





PREFACE 


The studies described in this report were performed principally by 
G. A. Marxman (theory), C. E. Wooldridge (theory and experiments), and 
R. J. Kier (experiments). The overall program has been under the direc- 
tion of G. A. Marxman. 

Administration and technical direction of the program have been under 
W. P. Peck, High Temperature Materials Branch, Langley Research Center. 

The authors wish to acknowledge the major contributions of A. J. Amaro , 
R. G. McKee, Jr., and W. H. Johnson to the experimental program. The com- 
puter calculations described in this report are the work of Miss Margery 
Brothers . 


v 




CONTENTS 


ABSTRACT . , , , , , . iii 

PREFACE V 

LIST OF ILLUSTRATIONS ix 

NOMENCLATURE xiii 

I INTRODUCTION , , , 1 

1 1 BACKGROUND 5 

III THEORETICAL STUDIES 7 

IV EXPERIMENTAL STUDIES 11 

A. Measurements of Regression Rate Behavior in the Pressure- 

Sensitive Regime 11 

B. Experiments with Variable Mass Flow and Variable Chamber 

Pressure . 25 

C. Instability Observations ....... 29 

D. Regression Rate Measurements with Fluorine ....... 34 

E. Regression Rate Measurements with Nitrogen Tetroxide . . 41 

F. Experiments with an Orifice in the Grain 45 

G. Slab Burner Experiments ....... 49 

H. Experiments in the Combustion Simulator ......... 51 

V EXPERIMENTAL RESULTS AND COMPARISON WITH THEORY 55 

VI CONCLUSIONS AND RECOMMENDATIONS 63 

APPENDICES 

A DIFFUSION-LIMITED THEORY OF STEADY-STATE HYBRID COMBUSTION . 65 

B GAS-PHASE EMISSIVITY 73 

C DIFFUSION-LIMITED REGRESSION RATE CALCULATIONS 77 

D DEVELOPMENT OF A MODEL FOR THE PRESSURE-SENSITIVE DOMAIN . . 89 

vii 



CONTENTS 


E REGRESSION RATE MEASUREMENTS WITH OXYGEN 93 

F CHAMBER PRESSURE RESPONSE TO VARIABLE OXIDIZER MASS FLOW . . 99 

G THERMOCHEMICAL CALCULATIONS FOR HYBRID PROPELLANTS 101 

REFERENCES 115 


viii 



ILLUSTRATIONS 


1 Influence of Activation Energy on Regression Rate Behavior , 8 

2 Influence of Gas-Phase Order on Regression Rate Behavior . . 9 

3 Weight Loss of a Pure PBAN Hybrid Grain (D = 1 in,) .... 13 

o 

4 Weight Loss of a PBAN Hybrid Grain Containing 20 Percent 

Aluminum (D = 1 in.) . . 14 

o 

5 Weight Loss of a PBAN Hybrid Grain Containing 40 Percent 

Aluminum (D = 1 in.) 15 

o 

6 Weight Loss of a Pure PBAN Hybrid Grain (D =1.5 in.) ... 16 

o 

7 Weight Loss of a PBAN Hybrid Grain Containing 20 Percent 

Aluminum (D = 1.5 in.) 17 

o 

8 Weight Loss of a PBAN Hybrid Grain Containing 40 Percent 

Aluminum (D = 1.5 in.) ... 18 

o 

9 Typical Chamber Pressure Traces from the Firing of a Pure 

PBAN Hybrid Grain 19 

10 Typical Chamber Pressure Traces from the Firing of a PBAN 

Hybrid Grain Containing 40 Percent Aluminum .... 20 

11 Experimental Dependence of Regression Rate on Pressure for 

the Aluminized PBAN-Oxygen Hybrid System (D = 1 in.) . . , . 20 

o 

12 Experimental Dependence of Regression Rate on Pressure for 

the Aluminized PBAN-Oxygen Hybrid System (D =1.5 in.) ... 21 

o 

13 Weight Loss of a Pure PU Hybrid Grain (D = 1 in.) 22 

o 

14 Weight Loss of a PU Hybrid Grain Containing 20 Percent 

Aluminum (D = 1 in.) 23 

o 

15 Weight Loss of a PU Hybrid Grain Containing 40 Percent 

Aluminum (D = 1 in.) 24 

o 

16 Chamber Pressure Response of the Hybrid Test Motor to Oxi- 
dizer Flow Perturbations 26 

17 Layout of Flow System for Oscillatory Oxidizer Flow Studies . 28 

18 Drive Mechanism for Oscillatory Chamber Pressure Studies , , 28 

ix 



ILLUSTRATIONS 


19 Inherent Hybrid Instability Wave Observed at Low Pressure 

in PBAN-Oxygen Propellant System 30 

20 Effect of Mass Flow Variations on Low Pressure Instability in 

Pure PBAN-Oxygen Propellant System ........ 32 

21 Instability Wave Observed at High Pressure in Pure PBAN- 

Oxygen Propellant System , 33 

22 Fluorine Hybrid System Assembly ... 35 

23 Internal Assembly of Fluorine System .... 36 

24 Injector Used with the Fluorine System 37 

25 Schematic of Fluorine Flow System 38 

26 Weight Loss of a Pure PU Hybrid Grain Burned with Fluorine 

(D = 1 in.) 39 

o 

27 Weight Loss of an 80 Percent PU/20 Percent Aluminum Hybrid 

Grain Burned with Fluorine (D = 1 in.) 40 

o 

28 N O Hybrid System Assembly 41 

£ * 

29 Injector Used with the N^O^ Hybrid System .......... 42 

30 Water Spray Pattern from N^O^ Injector ..... 42 

31 Schematic of N O Flow System 43 

2 4 

32 Weight Loss of a Pure PU Hybrid Grain Burned with N O 

2 4 

(D = 2.8 in.) 44 

o 

33 Weight Loss of a Pure PU Grain Burned with NO ....... 45 

2 4 

34 An Orifice Installed in a PU Grain 46 

35 Weight Loss of a Pure PU Grain Containing an Orifice and 

Burned with Oxygen (D =1.5 in.) 47 

o 

36 Spontaneous Instability Observed with an Orifice in the Grain 

(Pure PU) 48 

37 Small-Scale Slab Burner 49 

38 Schlieren Photographs of a Burning Plexiglas Slab 50 


x 



ILLUSTRATIONS 


3 9 Hybrid Combustion Simulator . . 52 

40 Temperature Profiles in a Burning Turbulent Boundary Layer . 54 

41 Regression Rate Behavior of PU-Oxygen Propellant System ... 57 

42 Regression Rate Behavior of PBAN-Oxygen Propellant System . , 59 

43 Effect of Mass Flux on Regression Rate 61 

44 Regression Rate Behavior of PU-Fluorine Propellant System . . 61 

A-l Theoretical Hybrid Mass Transfer Number .... 72 

B-l Coordinate System for Calculation of Emissivity . . 73 

C-l Variation of Mass Transfer Number with Aluminum Loading for 

PU and PBAN Binders 77 

C-2 Summary of Input Required for Hybrid Computer Code ..... 79 

C-3 Hybrid Regression Rate Computer Code 82 

C-4 Typical Output Format for Computer Calculations of Hybrid 

Regression Rates . 87 

E— 1 Small-Scale Tubular Hybrid Motor ..... 94 

E-2 Oxygen Flow-Control System 95 

E-3 Weight Loss of a Plexiglas Grain as a Function of Chamber 

Pressure (D = 1 in.) 97 

o 


xi 




NOMENCLATURE 


Port area for flow 
Nozzle throat area 

Thermochemical mass transfer number 

2 (p v) /p u c 
w e e f 

„ 2 

Friction coefficient = 2t /p u 

w e e 

Constant 

Stanton number = Q /p u Ah 
w e e 

Heat capacity of nonvaporizing fuel component 
Specific heat at constant pressure 
Specific heat at constant volume 
Diameter 

Gas-phase activation energy 
Total mass flux per unit area 

Oxidizer mass flux per unit area at x = 0 and t 
Enthalpy 

Heat of gasification of solid phase 


Radiative intensity per unit volume 



Thermal conductivity of the gas phase 


Mass fraction of nonvaporizing material in the solid phase 

Mass fraction of volatile surface material which forms solid combus- 
tion products; chamber parameter 

Chamber parameter 

Characteristic flame propagation length 

Characteristic turbulent mixing length 

Length of grain 
Lewis number 

Heat of vaporization of solid phase 

Mass flow rate 
Molecular weight 

Particle number density per unit volume of gas; order of the gas- 
phase reaction 

Pressure 

Perimeter of grain 

Prandtl number 

Chemical heat of reaction 

Convective heat transfer in the absence of radiation 
Radiative heat transfer 
Heat transfer to the wall 


Regression rate 



r Radius of a radiating particle 
P 

R Outer radius of radiating annulus; gas constant 

R Inner radius of radiating annulus 

o 

Re Reynolds number based on x 

x 

t Time 

t Characteristic time associated with combustion 

T Temperature; characteristic time associated with turbulent mixing 
u X-component of velocity 

u^ Kinetic flame speed 

u 7 Characteristic turbulent fluctuation 
Up Flame temperature 

V Volume 

x Longitudinal coordinate in stream direction 
x A Longitudinal coordinate of boundary layer merge point 

X Particle displacement by turbulence 

y Coordinate normal to solid surface 

z Radiative path length 

a Empirical radiation coefficient 

g Angle defined in Fig. 2 

Y Specific heat ratio, c /c ; absorptivity coefficient 

P v 

6* Boundary layer displacement thickness 

xv 



Ah 


Sensible enthalpy difference between combustion zone and wall 


e Emissivity 

Q Mass of oxidizer consumed in producing particulate products per unit 
mass of nonvaporizing fuel particles 




Mass of oxidizer consumed in producing particulate combustion prod- 
ucts per unit mass of volatile surface material which forms particu- 
late products 


\ Mass of particles produced per unit mass of particle-producing ma- 
terial in surface 


v Kinematic viscosity 


v Stoichiometric coefficient 
P 


v Stoichiometric coefficient 
r 


p Density 

n Particle density 

K P 

<j Stef an-Bol tzmann constant 

T Shear stress; transmissivity 

9 Angle defined in Fig. 2 

$ Total radiation per unit time 
Subscripts 

b Condition at combustion zone 

e Conditions at boundary layer edge 

f Fuel 

g Gas phase 

o Conditions with zero surface mass injection; oxidizer 


xv i 



r Radiative 

v Conditions of vaporizing component in solid phase 

w Conditions at surface 

Limiting value at high pressure 


xvii 




INTRODUCTION 


I . 

Increased interest in the development of the hybrid rocket as an im- 
portant propulsion device has naturally led to a consideration of its 
behavior over a wide range of operating chamber pressures. Earlier studies 
of hybrid combustion have concentrated on the high-pressure regime in which 
the regression rate is controlled by the turbulent diffusion mechanism in 
the boundary layer formed over the vaporizing surface. For this operating 
regime a relatively comprehensive theoretical description has been devel- 
oped. 1-6 For most propellants the theoretical model is in good agreement 
with experimental observations at pressures above about 200 to 300 psi , 7 > B 
The theory shows that in diffusion-limited combustion the regression rate 
is independent of pressure for a given mass flux through the port. 6 > 8 

For each propellant there is a pressure threshold above which the 
diffusion-limited theory becomes valid. Below this threshold there is a 
decrease in the regression rate with a decrease in the pressure at a fixed 
oxidizer mass flux. 9 > 10 The next step in the development of hybrid com- 
bustion theory must be a description of the low-pressure regime, in which 
the regression rate is pressure sensitive. 

Smoot and Price attempted to correlate their data in a semiempirical 
fashion, assuming heterogeneous attack to be the rate-controlling factor. 11 
An analysis by one of the present authors, however, has shown the hetero- 
geneous reactions alone cannot lead to the observed pressure dependence. 8 

Miller has proposed a model in which the pressure dependence is at- 
tributed to the broadening of the gas-phase flame zone caused by reduced 
reaction rates. 13 Several objections may be raised to this treatment. 

For example, isothermal conditions are assumed , and convective transfer 
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is neglected relative to diffusive transport. These assumptions are not 
applicable to hybrid combustion. 

A third treatment, advanced by Kosdon and Williams, is an analysis 
of small departures from the pressure-independent, diffusion-limited re- 
gime. 13 Comparisons of this theory with measurements taken under condi- 
tions far from the diffusion-limited regime, for which the theory cannot 
be valid, show much better agreement than in the small-departure regime, 
to which the theory should be applicable. 13 

The presence of pressure sensitivity in the regression rate leads to 
the possibility of an acoustic combustion instability mechanism similar 
to that encountered in solid propellants. A review of the existing liter- 
ature shows that little of the research conducted in this country has been 
directed specifically toward the combustion instability problem in hybrid 
motors. A study of available instability data suggests, however, the ex- 
istence of at least three major categories of hybrid combustion instability 
phenomena : 

1. Instability due to the lag time associated with vaporization and 
combustion of liquid droplets (with liquid injection) . This type 
of instability is related to that in liquid propellant motors but 
is probably considerably more complex, owing to the influence of 
the solid-liquid interface processes that affect lag times. 

2. Unstable burning associated with periodic accumulation and break- 
off of char layers or melted layers at the surface. The result- 
ing instability is generally low in frequency (like "chuffing") 
and occurs during operation at low regression rates. 

3. The combustion instability mechanism mentioned previously, which 
is related to that encountered in solid propellants. Here the 
acoustic admittance of the reacting turbulent boundary layer plays 
a key role, and the instability becomes more severe with a pres- 
sure-sensitive regression rate. 

The first year's effort on this program has been devoted largely to the 
third type of instability, which is potentially the most critical. Because 
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the presence of acoustic instability depends on the presence of a pressure- 
sensitive regression rate behavior, an understanding of the steady-state 
combustion mechanism that leads to pressure sensitivity in this regime is 
a prerequisite to development of a combustion instability theory. For 
this reason the main objective of the first year's effort has been to de- 
fine, both experimentally and theoretically, the dependence of regression 
rate on pressure in the kinetically controlled regime for several repre- 
sentative hybrid propellant combinations. It was found that a relatively 
simple analytical model based on classical turbulent flame theory exhibits 
good agreement with the observed regression rate/pressure dependence. 
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II. BACKGROUND 


Before proceeding to a discussion of the work carried out under the 
current program, it is worthwhile to review the previous status of hybrid 
combustion fundamentals. The basic premise underlying existing theory 
is that the combustion process may be described through a boundary layer 
model. It has been well established 1 ' 3 that in a boundary layer with com- 
bustion and a mass addition rate typical of the hybrid application, trans- 

4 

ition occurs at a length Reynolds number, Re of approximately 10 , or very 
near the leading edge of the grain. Therefore, a turbulent boundary layer 
treatment is required for a heat transfer theory of hybrid combustion. 

All present treatments of turbulent flow are necessarily semiempirical , 
but the cumulative efforts of a number of investigators such as Lees 14 and 
Spalding 16 have resulted in a reasonably satisfactory and remarkably simple 
description of heat and mass transfer in turbulent boundary layers with 
chemical reactions. 

A detailed development of the theory describing hybrid combustion in 
the pressure-insensitive (or high-pressure) regime is presented in Ap- 
pendix A. It is shown that in a hybrid motor in which radiative heat trans- 
fer is negligible, the regression rate is given by: 

Q 

c - 0.6 .0.23 -0.2 

p r = - — = 0.036G (p/p ) b' Re (1) 

f h ^ e x 

v 

The radiative heat transfer to the solid grain is given by: 

Q = cj e e (T 4 - T 4 ) (2) 

r w g r w 


5 



When radiation is present to a significant degree, as in a motor with 
high metal loading, the regression rate is described by the expression: 


P 



(3) 


where Q /h is given by Eq . 1 . It will be noted that when radiation is 
c v 

low i e,, when Q /Q « 1 , Eq . 3 reduces to Eq . 1. Equation 3 also shows 
’ r c 

the trade-off between convective and radiative transfer which occurs be- 
cause of the decrease in convective heat transfer with increasing surface 
mass addition. Procedures required to use these equations, including the 
method of evaluating the effective heat of gasification, h , and the mass 
transfer number, B, for a given propellant, are discussed in Appendix A. 
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III. THEORETICAL STUDIES 


Theoretical studies in this program were divided into three parts: 

(1) derivation of an expression for the gas-phase emissivity e which ap- 

g 

pears in Eq . 2; (2) calculation of asymptotic regression rates at high 
pressure, using Eq . 3, for the propellant systems chosen for experimental 
study during this investigation; and (3) formulation of a theoretical de- 
scription of regression rate behavior in the low-pressure, kinetics-limited 
domain. Parts 1 and 2 were required to complete part 3, because the theory 
of the latter expresses the regression rate in terms of the asymptotic 
limit obtained in part 2. Details of the studies are presented in Appen- 
dices B, C, and D. 

The main theoretical result of the first year's investigation was a 
combustion model describing the regression rate behavior in the pressure- 
sensitive regime. From this model the following equations were derived 
(Appendix D) : 



J i n/2 1 + n/2 - E /2RU 

2 P U e f f 


(4) 


(5) 


The constant C is related to flame speed kinetics. The characteristic 

length SL is associated with the scale of the combustion zone and JL is 
1 2 
identified with the scale of turbulence. At present the coefficient 

must be treated as an empirically determined constant for a given 
propellant system. 
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Typical results obtained from Eq . 4 are shown in Figs. 1 and 2. These 

results have been obtained by setting 0 = 1 at 10 atm pressure (to evaluate 

the parameter Ci / SL ) and using the thermochemical flame temperature versus 
1 2 

pressure behavior of a hybrid propellant system consisting of a 60 percent 
polyurethane/40 percent aluminum fuel with oxygen as the oxidizer. Fig- 
ure 1 shows the effect of changes in the gas-phase activation energy on 
the regression rate behavior when the order of the gas-phase reaction is 
held constant. Figure 2 shows the increased effect of changes in n with 

E held constant. It can be seen that n is a more powerful parameter than 
f 

E in affecting the shape of the regression rate curve. The proper choices 
for these parameters, as well as for C i/i must be found by comparison 
with experimental data. 

When the validity of Eq . 4 has been established by comparison with 
experimental data, it will be incorporated into the hybrid regression rate 
computer code for use during the second phase of the program. 



FIG. 1 INFLUENCE OF ACTIVATION ENERGY ON REGRESSION 
RATE BEHAVIOR 
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FIG. 2 INFLUENCE OF GAS-PHASE ORDER ON REGRESSION 
RATE BEHAVIOR 
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IV. EXPERIMENTAL STUDIES 


The experimental effort during the first year of the program was di- 
rected toward three objectives: (1) to delineate in detail the regression 

rate behavior in the pressure-sensitive regime; (2) to begin defining the 
combustion instability mechanism in the same regime; and (3) to begin de- 
veloping basic data on the structure of the boundary layer combustion proc- 
ess for eventual incorporation into an instability theory. During this 
year most of the effort centered on the first objective. This was the 
essential first step toward understanding combustion instability; before 
the latter can be described theoretically, it is necessary to understand 
the mechanism underlying regression rate dependence on pressure in the low- 
pressure regime, since instability is likely to result from the same mech- 
anism. Studies on objectives 2 and 3 were relatively preliminary in char- 
acter and will receive greater emphasis in future work, when the nature of 
combustion instability is investigated in more detail. 

A, Measurements of Regression Rate Behavior in the Pressure-Sensitive 
Regime 

Polyurethane and polybutadiene-acrylic-nitrile binders with varying 
aluminum loading were chosen as fuels and oxygen was used as the oxidizer 
in most of these tests. PU and PBAN were chosen for comparative binder 
tests because their surface behavior during regression is quite different. 
PU sublimes cleanly at the regressing surface, while the regression of PBAN 
consists of the stripping of hydrogen from the surface, leaving a carbona- 
ceous char layer that subsequently breaks away in particular form. 

Tubular burners of two sizes were used during the experimental stud- 
ies. The smaller of the two was designed to accept 2-it-in . O.D. grains 
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whose internal port size and length could be chosen to suit the require- 
ments of any given test. Two initial port sizes, 1 in. and 1-j in. in di- 
ameter, were used. This burner was originally designed to operate with 
gaseous oxidizers (normally oxygen) and incorporated a tubular entrance 
section 12-in. long which is designed to deliver a smooth, uniform oxidizer 
flow to the head end of the grain . It can be operated at pressures up to 
600 psia. During the later phases of this part of the investigation, this 
burner was modified for use with liquid fluorine. 

A larger motor case that accepted 5-in. 0.D, grains was also fabri- 
cated. This burner incorporated a liquid oxidizer injector, but may also 
be used with gaseous oxidizers if desired. Nitrogen tetroxide (N^O^) was 
used as the oxidizer for the tests carried out with this motor. It can 
also be operated at pressures up to 600 psia. 

Appendix E provides a detailed description of the apparatus and pro- 
cedures used in these experimental studies. 

Initial data from the PU and PBAN propellant systems were concentrated 

2 

at an oxidizer mass flux of approximately 0.1 lb/in. -sec. Data were ob- 
tained from 1-in.- and 1^-in .-diameter motors at four pressure levels rang- 
ing from about 30 psia to about 260 psia (the exact pressure level for each 
firing is given in the figures that follow) . Even though the nozzle was 
not always choked at the lowest pressure (the critical pressure ratio 
across the nozzle is 1.77 for y = 1.2 corresponding to a chamber pressure 
of 26 psia), measurement of the regression rate is significant. 

Figures 3 through 5 show the weight loss behavior as a function of 
time for the PBAN binder system containing 0, 20, and 40 percent aluminum 
in a 1-in. I.D. configuration. Similar data obtained from grains having 
an initial I.D. of 1-^ in. are shown in Figs. 6 through 8. In each case 
theoretical calculations are shown for the highest pressure level. 
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0 2 4 6 8 10 12 14 16 


TIME — sec 


TA-6638-6B 


FIG. 3 WEIGHT LOSS OF A PURE PBAN HYBRID GRAIN (D Q = 1 in.) 


It can be seen that there is a significant drop in weight loss (and 
therefore regression rate) as the pressure is lowered at all three metal 
loadings. At the highest pressure the weight loss agrees well with the 
calculated value based upon diffusion-limited theory. 

Typical chamber pressure traces measured at the head end of the grain 
for two of the PBAN systems at high and low pressure are shown in Figs. 9 
and 10. It can be seen that the pressure trace from the aluminized grain 
is more ragged at high pressure. The initial small pressure rise noted 
during the first stages of the high-pressure runs is caused by the preflow 
during the ignition period. 

A more interesting presentation of the regression rate data can be 
obtained by plotting the weight loss at a given pressure divided by 
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WEIGHT LOSS 







WEIGHT LOSS 







GB-6638-5 


FIG. 6 WEIGHT LOSS OF A PURE PBAN HYBRID GRAIN (D Q = 1.5 in.) 

the weight loss at the highest pressure as a function of pressure. At any 
given time this curve will depict the behavior of the average regression 
rate up to that time. Results for the 1-in,- and 1-|— in .-diameter PBAN 
runs are plotted in this manner in Figs. 11 and 12. The data show a sim- 
ilar pressure dependence for both diameters and show that the regression 
rate decays to approximately 0.4 of its high pressure value at 1 atm. 

It is also evident from Figs. 11 and 12 that there is a scatter of 
approximately ±10% in the data points obtained from weight loss measurements. 
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WEIGHT LOSS 



GB-6638-6 


FIG. 7 WEIGHT LOSS OF A PBAN HYBRID GRAIN CONTAINING 20 PERCENT 
ALUMINUM (D q = 1.5 in.) 
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TIME — sec 

GB-6638-7 

FIG. 8 WEIGHT LOSS OF A PBAN HYBRID GRAIN CONTAINING 40 PERCENT 
ALUMINUM <D q = 1.5 in.) 
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FIG. 9 TYPICAL CHAMBER PRESSURE TRACES FROM THE FIRING OF A PURE 
PBAN HYBRID GRAIN 

The agreement between theory and experiment obtained at high pressure 
shown in the previous figures is much closer, indicating that most of the 
scatter makes its appearance in the kinetically dominated regime. In 
general, however, regression rate data obtained from weight loss measure- 
ments exhibit scatter because they represent the quotient of two measured 
quantities, with the inherent possibility of doubling the experimental 
errors , 

Finally, in Figs, 11 and 12, both of the regression rates averaged 
over 5 and 10 sec exhibit nearly identical behavior, as would be expected 
since the change in oxidizer mass flux between the two times is not sig- 
nificant , 
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FIG. 10 TYPICAL CHAMBER PRESSURE TRACES FROM THE FIRING OF A PBAN 
HYBRID GRAIN CONTAINING 40 PERCENT ALUMINUM 
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FIG. 11 EXPERIMENTAL DEPENDENCE OF REGRESSION RATE ON PRESSURE FOR THE 
ALUMINIZED PBAN-OXYGEN HYBRID SYSTEM (D 0 = 1 in.) 
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64 - 6638-11 

FIG. 12 EXPERIMENTAL DEPENDENCE OF REGRESSION RATE ON PRESSURE FOR THE 
ALUMINIZED PBAN-OXYGEN HYBRID SYSTEM (D q = 1.5 in.) 

Weight loss data for aluminized PU grains are given in Figs. 13 
through 15. These data show a similar dependence of regression rate on 
pressure to that observed in the PBAN system discussed previously. For 
this reason no tests were run in a l-|-in . initial port diameter configura- 
tion, Results of tests with varying flow, which are discussed in detail 
below, are also shown in Figs. 13 through 15. The dependence of the re- 
gression rate of the PU grains on pressure will be discussed in a later 
section . 
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FIG. 13 WEIGHT LOSS OF A PURE PU HYBRID GRAIN (D Q = 1 in.) 
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B. Experiments with Variable Mass Flow and Variable Chamber Pressure 

In the pressure-sensitive domain of hybrid combustion the pressure 
sensitivity may give rise to an acoustical instability, characterized by 
the "acoustical admittance," which is defined as the ratio of the regres- 
sion rate perturbation induced by a given pressure perturbation to the 
pressure perturbation, or Ar/Ap. As pointed out earlier, the derivation 
of an analytical expression for the acoustic admittance depends on the 
initial formulation of an expression for the steady-state regression rate/ 
pressure dependence; this has been the object of the theoretical and ex- 
perimental work described above. 

Experiments to be described below were conducted to define the con- 
ditions under which hybrid instability arises and to delineate the coupling 
mechanism. It is well known that solid propellants exhibit unstable be- 
havior at frequencies of the order of hundreds of cycles per second. 16 , 19 
Since the regression rates of the hybrid systems being considered here 
are only on the order of one-quarter of typical solid propellant burning 
rates, and since the characteristic frequency of maximum response depends 
on the square of the burning rate, it was decided to conduct a few pre- 
liminary oscillatory experiments by oscillating the oxidizer flow and the 
nozzle throat area at about 10 cps while observing the resulting chamber 
pressure. This procedure makes it possible to distinguish different coup- 
ling mechanisms. 

The oxidizer mass flow can be held constant during nozzle throat area 
variations by employing a sonic choke at the inlet. On the other hand, 
mass flux perturbations that are introduced by oxidizer flow variations 
will give rise to pressure perturbations that are induced by the varying 
mass flow through the nozzle. In fact, these pressure perturbations can 
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be used as a measure of the mass flux perturbations. Obviously, great care 
must be taken in the subsequent interpretation of results to carefully 
separate pressure and mass flux effects. 

Another facet which must be carefully considered is the chamber re- 
sponse. An oscillating oxidizer mass flow will not introduce a correspond- 
ing mass flux oscillation in the chamber if the chamber response function 
is highly damped at the frequency of interest. Obviously, the chamber 
response function must be understood and known before the regression rate 
response can be correctly deduced. This point is discussed in detail in 
Appendix F. In Fig. 16 the calculated chamber response is shown for 
two cases, isothermal and isentropic behavior. The isothermal response 
is of course somewhat lower than the isentropic response, owing to losses 
in the isothermal case. 



NON DIMENSIONAL FREOUENCY, rw 

FIG. 16 CHAMBER PRESSURE RESPONSE OF THE HYBRID TEST MOTOR 
TO OXIDIZER FLOW PERTURBATIONS 
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The salient fact that emerges from a consideration of Fig. 16 is that 
the pressure response is very low at the higher values of frequency. Since 
this response is a measure of the mass flux perturbation in the chamber, 
it is evident that at, say, 35 cps a very large upstream oxidizer mass 
flow perturbation will be required to obtain a relatively small mass flux 
perturbation in the chamber. 

A similar calculation could be made for the variable pressure- 
constant mass flux case. The general shape of the solution would show 
a similar frequency dependence. It is clear that care must be taken to 
avoid attributing chamber response effects to the propellant response 
which is of main interest. 

The experimental systems used for the oscillatory oxidizer flow and 
chamber pressure studies are shown in Figs. 17 and 18. A ball valve which 
is driven by a variable speed motor was inserted in the oxidizer flow line 
just above the inlet sonic choke (see Fig. 17) to provide flow variations. 
The frequency of oscillation is determined by the speed of the motor and 
the peak amplitude by the upstream pressure. To achieve chamber pressure 
oscillations separately, a nozzle slug assembly was built as shown in 
Fig. 18. The maximum frequency obtainable from both the ball valve and 
the oscillating probe was limited to about 30 cps by mechanical inertia. 

Two experimental results are of particular interest in these two 
tests: first the determination of whether the mean regression rate (cor- 
responding to the mean pressure) is changed by the presence of the oscil- 
lation, and second, the comparison of the observed magnitude of the chamber 
pressure oscillation to the computed magnitude, which is based upon a 
quasi-steady relationship between regression rate and pressure. 
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FIG. 17 LAYOUT OF FLOW SYSTEM FOR OSCILLATORY OXIDIZER 
FLOW STUDIES 



FIG. 18 DRIVE MECHANISM FOR OSCILLATORY CHAMBER 
PRESSURE STUDIES 


28 






Weight loss data obtained in the PU system under conditions of varying 
flow were previously shown in Figs. 13 through 15. The mass flux was 
varied approximately by a factor of 2 at a frequency of about 10 cps . It 
can be seen that the variable flow rate data follow the same trends as 
the constant flow rate data with no discernable differences. Similar re- 
sults were obtained in the PBAN system and for pressure changes induced 
by a variable throat area. A complete summary of all mean regression 
rate data is presented in a later section. 

For a typical run in which the mass flux was varied between 0.075 and 
2 

0.15 lb/in. -sec (with an average value of 0.11) at 12 cps, the chamber 
pressure was observed to oscillate with an amplitude of ± 7 psi at an 
average value of 43 psia. Using representative values of V = 20 in . and 
T = 5000°R, the calculated value of the nondimensional frequency tul) is 
5.2 and the value of the response is 0.35. As shown in Fig. 16, this re- 
sponse is approximately 50 percent higher than the calculated quasi-steady 
value. Thus, although the mean regression rate is undisturbed by the oscil- 
lating mass flux, evidence of a propellant response is apparent. Further 
interpretation of this data must await the theoretical analysis of the 
transient regression rate behavior, which is planned for the second year 
of the program. 

C . Instability Observations 

The mean regression rate data reported above, obtained for both 1-in.- 
and 1 . 5-in , -diameter grains, exhibited a similar mean regression rate/ 
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pressure dependence for both grain sizes. However, during the burning of 
the 1 .5-in .-diameter grains, an instability wave was observed to develop 
when the tube diameter reached approximately 1.75 in. This event appears 
to be an inherent instability in the combustion process. 

Typical wave forms for the PBAN-oxygen propellant system operating 
at a mean pressure of 2.5 atm are shown in Fig. 19. The upper trace in 



(c) 80% PBAN 

20% ALUMINUM 


p = 2.5 atm 


(d) 60% PBAN 

40% ALUMINUM 

TA -6638 20 


FIG. 19 INHERENT HYBRID INSTABILITY WAVE OBSERVED AT LOW PRESSURE 
IN PBAN-OXYGEN PROPELLANT SYSTEM 
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each chart is the energy in the frequency band between 150 and 1500 cps . 
At this mean pressure level the mean regression rate is approximately 0.6 
of the asymptotic value at high pressure. The frequency of the predomin- 
ant mode is about 650 cps, and the peak-to-peak amplitude of the wave is 
of the order of 10 psi , or 25 percent of the mean pressure. This wave 
form has been observed with all metal loadings up to 40 percent, as shown 
in Fig. 19. Such behavior is consistent with the previous observation 
that all formulations exhibited the same regression rate/pressure depend- 
ence. It is not surprising that the presence of metal particles has no 
damping effect, because the observed frequency is too low to be affected 
by particle damping phenomena. Of perhaps more fundamental importance is 
the observation that the addition of burning metal particles does not in- 
fluence either the dependence of the mean burning rate on pressure or the 
frequency of the observed instability wave. Such behavior implies that 
the burning process (i.e., flame speed) of the volatile products released 
by the vaporization of the binder is of primary importance in establishing 
the pressure-dependent regime of hybrid combustion. 

The rather low frequency which is observed indicates the excitation 
of a longitudinal rather than a transverse mode. In this connection it 
was found during the studies with variable mass flow that a very low fre- 
quency variable mass flow had a definite effect on the instability process. 
Typical results for the unfilled PBAN binder system are shown in the se- 
quence of photographs presented in Fig. 20. The mass flux was varied 
30 percent at a frequency of approximately 10 cps. Photograph (a) shows 
that on the average the disturbance is magnified near peak values of the 
mass flux (or chamber pressure, since the two are interrelated) and damped 
near trough values. Photographs (b) and (c) show the behavior near a 
peak at progressively faster sweep rates. The upper trace is again fil- 
tered to show the energy between 150 and 1500 cps. In this case the 
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(a) 50 msec/cm SWEEP RATE 



(b) 10 msec/cm SWEEP RATE 



p = 6 atm 


(c) 2 msec/cm SWEEP RATE 

FIG. 20 EFFECT OF MASS FLOW VARIATIONS ON LOW PRESSURE 

INSTABILITY IN PURE PBAN-OXYGEN PROPELLANT SYSTEM 


instability amplitude is about 15 psi at a mean pressure level of about 
90 psia. The frequency is about 650 cps as before. 

Figure 21 shows the behavior of the unfilled PBAN binder system 
operating at 15 atm with a steady oxidizer flow rate. At this pressure 
level the mean burning rate is accurately predicted by the heat-transfer- 
limited theory, and a severe pressure-coupled instability would not be 
expected. Figure 21 shows that a small pressure disturbance of about 
10 psi (compared to a mean level of 225 psia) is present; this distui bance 
does not exhibit the pure waveform of the instability observed at lower 
pressures. It is also of interest to note that the lower unfiltered tiace 
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FIG. 21 INSTABILITY WAVE OBSERVED AT HIGH PRESSURE 
IN PURE PBAN-OXYGEN PROPELLANT SYSTEM 


contains a 2-psi wave at a frequency of about 3500 cps . This high-frequency 
wave may represent a transverse mode of instability. In all cases the in- 
stability arose spontaneously without external excitation. 

According to the thermochemical calculations presented in Appendix G 
the value of the specific heat ratio y is approximately 1.2 for all of the 
PU and PBAN propellants. The critical pressure ratio across the exit noz- 
zle required to maintain sonic flow at the throat is given by 


P 

c 

P 

a 



( 6 ) 


For y = 1.2, the value of this ratio is 1.77, giving p = 26 psia for 

c 

standard conditions. The lowest experimental mean pressure ratio was 2.5 
with a sinusoidal variation of ±0,3, Therefore, the observed behavior 
cannot be attributed to an unsteady choking and unchoking of the nozzle . 

Assuming that the wave is longitudinal in nature (i.e., an organ 
pipe oscillation) , the average wave speed can be calculated as twice the 
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length of the chamber (4 ft) times the measured frequency of the wave: 

u = 4(650) = 2600 ft/sec 
w 

At first glance, this speed appears to be lower than the speed of sound 
one would expect in the hybrid motor. However, two points must be borne 
in mind: First, the motor configuration used in these studies consisted 

of a 12-in. -long straightening section upstream of the grain. In this 
section the temperature, and therefore the speed of sound, is low. Second, 
the cool core in the motor itself will tend to decrease the average sonic 
speed. The result is that the observed wave speed is a complicated product 
of the environment. Future tests involving different chamber lengths will 
be carried out to determine the influence of the coupling between chamber 
length and propellant response. 

The question of fundamental interest is how the propellant response 
acts as a driving mechanism to support the wave under certain circumstances, 
i ,e . , with a certain motor diameter at a given oxidizer mass flux. This 
question will be investigated in detail during the future phases of the 
program . 

D. Regression Rate Measurements with Fluorine 

In addition to the many experiments described previously with oxygen 
as the oxidizer component of the propellant system, a few runs were made 
with liquid fluorine as the principal oxidizer with only enough oxygen to 
burn the carbon in the binder to CO. Two fuel systems were used, 100 per- 
cent PU and 80 percent PU/20 percent A1 . The objective of these experi- 
ments was to compare the low-pressure regression rate behavior of the 
fluorine propellant system with that of the oxygen propellant system. 

The total number of fluorine runs was limited by the operating costs of 
the cryogenic system. 
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A photograph of the test apparatus used in the fluorine runs is in- 
cluded as Fig. 22. The dewar flask contains a liquid nitrogen bath in 
which all components of the fluorine system up to the injector face are 
contained. Figure 23 shows the internal assembly, including the fluorine 
tank that is pressurized by helium and the coil that is used to liquify 
the fluorine. The injector consists of a simple central hole (0.055 in. 
in diameter) for the fluorine with a coaxial annulus (0.085 in. clearance) 
for the oxygen. Figure 24 is a diagram of the injector. The entire as- 
sembly up to but not including the pneumatic actuators for the valves is 
immersed in liquid nitrogen. Figure 25 is a schematic diagram of the 
fluorine system. The liquid fluorine, pressurized with gaseous helium 
is metered through a cavitating venturi. The injector is prechilled with 



FIG. 22 FLUORINE HYBRID SYSTEM ASSEMBLY 
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FIG. 23 INTERNAL ASSEMBLY OF FLUORINE SYSTEM 


liquid nitrogen and the nozzle is cooled with circulating water. Note 
that no flow-straightening section such as that used in the tests with 
gaseous oxygen is required upstream of the grain because the rapid evapora- 
tion and expansion of the fluorine fills the inlet section uniformly. 

Weight loss data obtained from aluminized PU grains (0 and 20 percent 
Al) , are given in Figs. 26 and 27. The ratio of fluorine to oxygen was 
approximately 3:1; only enough oxygen was used to convert the carbon in 
the binder to carbon monoxide. Runs were made for 3 and 6 seconds in- 
stead of the 5 and 10 seconds used in the oxygen case to limit the total 
heat transfer back to the cryogenic system. Also shown in the figures 
is the transport-limited regression rate calculated with the hybrid 
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FIG. 24 INJECTOR USED WITH THE FLUORINE SYSTEM 


computer code. It can be seen that the data at 50 psia and above agree 
reasonably well with the calculated regression rates, while the data in 
the vicinity of 20 psia tend to fall below the calculated values. 

The 50 psia data point in the 100 percent PU system and the 20 psia 
data point in the 80 percent PU/20 percent A1 system seem to be higher 
than the remaining points; unfortunately time was not available to repeat 
these runs but additional data will be obtained during future phases of 
the program. In general, the departure of the regression rate from the 
transport-limited value begins at a lower pressure with fluorine than 
with oxygen; a direct comparison is made in a later section. 
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FIG. 25 SCHEMATIC OF FLUORINE FLOW SYSTEM 
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FIG. 26 WEIGHT LOSS OF A PURE PU HYBRID GRAIN BURNED 
WITH FLUORINE (D. = 1 in.) 
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FIG. 27 WEIGHT LOSS OF AN 80 PERCENT PU/20 PERCENT ALUMINUM 
HYBRID GRAIN BURNED WITH FLUORINE (D n = 1 in.) 



E. Regression Rate Measurements with Nitrogen Tetroxide 


The larger 5-in .-diameter motor configuration has been assembled and 
a limited amount of initial data has been accumulated with nitrogen tetrox- 
ide (NO) used as the oxidizing species with 100 percent PU grains. A 
2 4 

photograph of the motor assembly is given as Fig. 28. Ignition is accom- 
plished with the injection of hydrozene (N^H^) during the first second of 


operation. The vortex injector that was used is shown in detail in Fig. 29. 
N^O^ is injected in a vortical pattern through eight holes of 0.0595 in. 
in diameter; the swirling flow impinges on the wall of the injector and 
breaks up into a fine mist. Figure 30 is a photograph of a typical water 
spray pattern. It is essential to obtain good mechanical breakup of the 
spray pattern in a noncryogenic system such as this; the vortex injector 



FIG. 28 N 2 0 4 HYBRID SYSTEM ASSEMBLY 
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FIG. 29 INJECTOR USED WITH THE N 2 0 4 HYBRID SYSTEM 


FIG. 30 



WATER SPRAY PATTERN FROM N 2 0 4 INJECTOR 
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works very well from this point of view. A schematic of the complete 
flow control system is presented as Fig. 31. 


During the first year of the program the objective of the NO studies 

£ * 

was merely to design and check out the system. It will be very useful 
in the study of low-frequency instabilities associated with droplet vapor- 
ization lags. The larger size of the motor also makes it valuable for 



o-X REMOTE VALVE 
REGULATOR 



FIG. 31 SCHEMATIC OF N 2 0 4 FLOW SYSTEM 
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scaling tests. Data obtained from two successive runs on one motor are 
shown in Figs. 32 and 33. The lower pressure was 100 psia (6.8 atm) at 
which the regression rate had fallen to about 0.8 of its asymptotic high- 
pressure value when oxygen was used as the oxidizer. In this case the 
decrease appears to be less than that, but more data are required to es- 
tablish the actual dependence. 



FIG. 32 WEIGHT LOSS OF A PURE PU GRAIN BURNED WITH N 2 0 4 
(D q = 2.8 in.) 


44 





FIG. 33 WEIGHT LOSS OF A PURE PU GRAIN BURNED WITH N„0 4 
(D c = 2.9 in.) 


F. Experiments with an Orifice in the Grain 

To gain additional information about the spontaneous instability de- 
scribed earlier, several PU grains were fired with graphite orifices cast 
in place one-quarter of the way down the grain from the head end. The 
initial internal diameter of the grain was 1.5 in. and the orifice diameter 
chosen was 0.75 in., giving a constriction ratio of 4. Figure 34 is a 
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FIG. 34 AN ORIFICE INSTALLED IN A PU GRAIN 

photograph of an orifice in place. Oxygen was used as the oxidizer in 
these tests . 

Two results were of interest in these tests — the effect of the orifice 
on the weight loss (or regression rate) and its effect on the observed in- 
stability wave. Weight loss data are shown in Pig. 35. It can be seen 
that the overall effect of the orifice is to increase the burning rate by 
at least 50 percent. The burning rate/pressure dependence, however, is 
similar to that observed in the system without the orifice, indicating that 
the effect of the orifice is one of gross turbulent mixing and is not a 
direct influence on the chemical kinetics. 

The spontaneous appearance of an instability wave, similar to that 
illustrated in Figs. 25 and 26, was observed even with the orifice in 
place. The results are shown in Fig. 36. The frequency of the wave is 
approximately 510 cps, compared to the 650 cps recorded earlier without 
the orifice. The amplitude of the wave is somewhat larger at the lowest 
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FIG. 35 WEIGHT LOSS OF A PURE PU GRAIN CONTAINING AN ORIFICE 
AND BURNED WITH OXYGEN (D Q = 1.5 in.) 

mean pressure, but is smaller at the highest pressure, compared to the re- 
sults obtained without an orifice. Other investigators 3 have used orifices 
to suppress the low-frequency instability associated with droplet vaporiza- 
tion lags; it is apparent that this device has little effect on the present 
acoustic mode of instability. Further investigation will be carried out 
as the program continues. 
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FIG. 36 SPONTANEOUS INSTABILITY OBSERVED 
WITH AN ORIFICE IN THE GRAIN 
(Pure PU) 
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G. 


Slab Burner Experiments 


A small-scale slab burner, depicted in Fig. 37, has been constructed 
and used to obtain initial schlieren photographs of the hybrid combustion 
process as an aid to the parallel theoretical effort . The inlet to the 
burner contains a flow-straightening section to assure smooth flow condi- 
tions in the test section . Glass windows are provided in the 2 in . by 2 in . test 
section. Two schlieren photographs of the combustion zone over a Plexiglas 
slab obtained at two pressure levels are included as Fig. 38. The flow is 
from right to left, and the schlieren light has been blocked off at the 

left edge of the picture to show the luminosity of the flame zone. For 

5 

these tests the Reynolds number was approximately 1.5 x 10 /in. Note that 
the position of the flame zone, as determined from both the luminosity and 
the change in the density derivative observed on the schlieren picture, 
is approximately 10 to 15 percent of the total boundary layer thickness 
away from the wall. This result agrees with previous measurements obtained 



FIG. 37 SMALL-SCALE SLAB BURNER 
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FIG. 38 SCHLIEREN PHOTOGRAPHS OF A BURNING PLEXIGLAS SLAB 


at atmospheric pressure. 1 Note also that the flame zone is much thicker 

at low pressure where kinetics play an important role. The two pictures 

o 

shown here were obtained at an oxidizer mass flux of 0.1 lb/in. -sec. 

Future studies will consider other propellants, other oxidizer mass fluxes, 
and additional levels of pressure. The observed position and thickness of 
the flame zone will be used to aid in the theoretical interpretation of 
combustion in the low-pressure regime . 


50 




H. Experiments in the Combustion Simulator 


An ultimate goal of this program is to develop a theoretical descrip- 
tion of the hybrid combustion instability mechanism, so that unstable burn- 
ing can be predicted to a reasonable degree. It is quite certain that this 
mechanism involves a complicated interaction between the turbulence struc- 
ture of the boundary layer and the reaction process. Therefore, to develop 
a realistic model it will be necessary to obtain some experimental data 
that help to reveal the real nature of the coupling mechanism. Past ex- 
perience has shown that it is virtually impossible to obtain such data in 
slab burners or laboratory motors. For reliable measurements a technique 
is required that accurately simulates hybrid combustion but provides con- 
trolled burning conditions and a fixed solid/gas interface for spatial 
reference. To fill this requirement, a wind tunnel having a test section 
7 in. by 7 in. in cross-section which contains a 6-in. -wide by 30-in. - 
long porous section in the upper wall has been built with Institute funds. 
The design incorporates honeycomb and screens in the settling chamber, 
followed by a 10:1 contraction ratio, to ensure smooth flow conditions in 
the free-stream of the test section. The similator, shown in Fig. 39, 
operates near atmospheric pressure. The size of the tunnel was dictated 
by the desire to obtain a thick enough boundary layer over the porous sec- 
tion to allow accurate probing of the layer. A mixture of hydrogen and 
nitrogen is injected as the fuel through the porous section for the com- 
bustion experiments. This fuel is used because of the wide flammability 
limits of the hydrogen-air flame which make it possible to vary the maximum 
temperature in the boundary layer over wide limits. The significance of 
this point will be discussed later. The boundary layer thickness which 
is obtained is approximately two inches, depending on the rate of wall 
mass injection. The test section has a moveable lower wall so that a wide 
range of axial pressure gradients can be imposed. The value of the simu- 
lator is that it allows detailed diagnostic measurements to be carried out 
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FIG. 39 HYBRID COMBUSTION SIMULATOR 


in a relatively thick, burning, turbulent boundary layer over a stationary 
(i .e . , nonregressi ve) surface. 

The theoretical analysis presented earlier requires as an input an 
estimate of the rms turbulence level in the presence of the flame zone. 

The functional dependence of the turbulence level on such parameters as 
mass flux, Reynolds number, and mass transfer number is of primary 
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importance in attaining a general understanding of the process. Current 
results, which are discussed in detail in the next section of this report, 
indicate that good agreement is obtained between the theoretical and ex- 
perimental regression rate/pressure dependence, but that the mass flux 
dependence is less satisfactory. 

Earlier turbulence measurements in a similar apparatus 20 showed a 
marked interaction between the flame zone and the turbulent structure, 
even though those measurements were restricted by the apparatus to a thin 
region near the wall . The indication to be gained is that further experi- 
ments closer to the flame zone with more sophisticated instrumentation will 
be required before any real understanding of the interaction is attained. 
For this reason, a constant-temperature hot wire anemometer, capable of 
being used with cooled probes, is being purchased for use during the second 
year of the program. 

Because hot wire measurements are required to define the fluctuation 
parameters of interest, only a limited number of experiments, consisting 
of temperature profile measurements, were carried out in the simulator 
during the first year's effort. Temperature profiles obtained in a zero- 
pressure gradient flow at two values of the mass transfer number B are 
shown in Fig. 40. 

Of particular interest is the observed change in the profile at a 
fixed value of B as the maximum temperature is reduced. This reduction 
is achieved by changing the proportion of hydrogen in the hydrogen-nitrogen 
fuel mixture injected through the porous plate. For the present experi- 
ments at B = 7,2 the reduction in maximum temperature from 950° C to 730° C 
was achieved by reducing the mass fraction of hydrogen in the incoming 
mixture from 4 to 3.5 percent. As the temperature is lowered the flame 
zone broadens considerably. From the point of view of the effect on 
turbulence structure in the boundary layer, such behavior simulates the 
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FIG. 40 TEMPERATURE PROFILES IN A BURNING TURBULENT BOUNDARY LAYER 

observed effects in the pressure-sensitive regime wherein the flame zone 
thickness grows as chemical kinetic processes begin to exert an influence. 
Turbulence measurements under the two conditions shown would help to es- 
tablish the dependence of the characteristic mixing time on the flame zone 
thickness and heat release distribution. 

Also shown in Fig. 40 is a profile obtained at B = 12.9. Note that 
this profile nearly matches the high-temperature profile for B = 7.2. 
Turbulence measurements with the higher rate of wall mass injection would 
also be valuable as an aid to the theoretical analysis of turbulent mixing 
and combustion. Such experiments are planned for the second phase of the 
program. 
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V. EXPERIMENTAL RESULTS AND COMPARISON WITH THEORY 


A complete synopsis of the experimental mean regression rate data is 
presented in Figs. 41 and 42 for the PU-oxygen and PBAN-oxygen propellant 
systems respectively. The results are presented as the mean regression 
rate over 5 sec of operation versus the mean pressure in logarithmic co- 
ordinates. The tests encompass three types: (1) constant oxidizer flow 

rate, constant nozzle area; (2) constant oxidizer flow rate, variable 
nozzle area; and (3) variable oxidizer flow rate, constant nozzle area. 

For the variable tests, the mean regression rate is plotted versus the 
mean pressure . 

First, it is observed from Figs. 41 and 42 that the presence of metal 
loading does not materially affect the observed pressure sensitivity when 

the results are presented in terms of the ratio r/r . (It should be noted 

00 

that r^ for the 40 percent aluminum-loaded systems is approximately twice 
that for the systems with no aluminum.) Also, within the scatter of the 
data, which seems to increase somewhat as the pressure is lowered, the 
variable mass flux and nozzle area tests (which both induced pressure vari- 
ations) gave results that are consistent with each other and with the re- 
sults from tests with constant conditions. 

Also shown in the figures are theoretical curves for the 100 percent 
binder and the 60 percent binder/40 percent aluminum fuel systems . These 
curves have been generated by assuming that t /T = 1 at 9 atm, n = 2, and 
E^ = 60 kcal/mole . The slight difference in the two curves at low pressure 
stems from the slightly different behavior of flame temperature with pres- 
sure for the two systems (see Appendix G) . For comparison, the least- 
square-error line through the data points is also shown. 
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FIG. 41 REGRESSION RATE BEHAVIOR OF PU-OXYGEN PROPELLANT SYSTEM 
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FIG. 42 REGRESSION RATE BEHAVIOR OF PBAN-OXYGEN PROPELLANT SYSTEM 



As shown in Section III, the exact choice of gas-phase activation 
energy is not critical to the calculated value of regression rate ratio, 
whereas the choice of gas-phase order n has a sizeable effect. The present 
choice of n = 2 gives very good agreement between the theory and the least- 
square-error curve through the data. This value of n implies that the pre- 
ponderance of heat release takes place in bimolecular reactions, as would 
be expected in a hydrogen/carbon/oxygen gas-phase flame. 

The data shown in Figs. 41 and 42 were obtained at an oxidizer mass 
flux of 0.1 lb/in. -sec. The data of Smoot and Price 9 ’ 10 show a similar 
dependence of regression rate on pressure at this mass flux and indicate 
further that the length ratio 4 /j& in Eq . 5 is a function of mass flux. 

This behavior is demonstrated by the data shown in Fig. 43. At high 

values of oxidizer mass flux and a fixed pressure, the analysis predicts 

0.4 

a regression rate proportional to G while the experimental regression 
rate becomes independent of the mass flux. It is not surprising that the 
unknown length ratio is some function of the mass flux because the trans- 
port of oxidizer through the boundary layer is a function of G. Further 
experiments at different values of mass flux must be carried out to settle 
this point. It is certainly true, however, that the theoretical prediction 
at a fixed mass flux agrees remarkably well with the data. 

The results from the test firings with fluorine are shown in Fig. 44. 
In this case the critical pressure, at which the parameter 0 is unity, is 
considerably lower than in the oxygen case, being about 2.25 atm instead 
of 9 atm (the theoretical curve for the oxygen system is shown for com- 
parison) . This behavior reflects the greater reactivity of fluorine com- 
pared to oxygen, thereby lowering the pressure threshold at which chemical 
kinetic processes begin to exert an influence. This implies that fluorine 
systems are less likely to exhibit acoustic instability; this point will 
be investigated as the program proceeds. 
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FIG. 43 EFFECT OF MASS FLUX ON REGRESSION RATE 



FIG. 44 REGRESSION RATE BEHAVIOR OF PU-FLUORINE PROPELLANT SYSTEM 
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VI. CONCLUSIONS AND RECOMMENDATIONS 


During the past year a major effort has been made to characterize the 
behavior, both experimentally and theoretically, of typical hybrid pro- 
pellant systems in the regime where the regression rate is pressure- 
dependent . The occurrence of a pressure dependence leads to the strong 
possibility of an acoustic instability, arising through the response of 
the regression rate to pressure perturbations. To obtain a mathematical 
model of the response function, it is first necessary to understand the 
steady-state pressure-sensitive behavior of hybrid propellants. 

Two binders were chosen for the experimental investigation — cleanly 
gasifying PU and charring PBAN to determine whether different surface 
kinetics play a first-order role in the combustion process. Aluminum 
loadings of 0, 20, and 40 percent were employed with both binder systems 
to ascertain the role of radiation in the pressure-sensitive regime. Oxy- 
gen, fluorine, and nitrogen tetroxide were used as oxidizers to help es- 
tablish the role of oxidizer reactivity in the combustion process. 

The main results of the first year's effort can be summarized as 
follows : 

1. The regression rate/pressure dependence observed at low pressures 
is independent of binder pyrolysis characteristics and aluminum 
loading for a given oxidizer, 

2. A noticeable decrease is observed in the regression rates of PU 
and PBAN fuels below 10 atm when oxygen is used as the oxidizer 
and below 2.5 atm when fluorine is used. These values pertain to 
an oxidizer mass flux of 0.1 lb/in. ^-sec and would decrease with 
increasing mass flux. The lower threshold pressure for fluorine 
is a consequence of its generally greater reactivity. 
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3 . A relatively simple analytical model based on the wrinkled-flame- 
sheet turbulent flame theory has been developed . The model agrees 
well with the data obtained at a fixed oxidizer mass flux. Modi- 
fications are needed to improve the description of variable mass 
flux conditions . 

4. A spontaneous longitudinal instability was observed to develop 
under certain conditions of diameter and mass flux. The presence 
of this wave is related to the ability of the regression rate to 
respond to a pressure perturbation of the given frequency. The 
wave was virtually unaffected by the addition of an orifice in 
the grain. 

5. Experiments with a pulsed oxidizer flow showed a greater varia- 
tion in the chamber pressure than could be calculated by assuming 
a quasi-steady regression rate. This result indicates a propel- 
lant response of the kind that arises in unstable combustion. 

The pulsed flow also acted to damp the higher frequency spontan- 
eous oscillation. 

6. Schlieren photographs obtained from the slab burner showed that 
the combustion zone in the turbulent boundary layer is much 
thicker at lower pressure where chemical kinetics exert a sig- 
nificant influence than at high pressure. 

During the second phase of the study, additional regression rate data 
will be obtained at other values of the oxidizer mass flux to provide scal- 
ing information for further theoretical development . Longitudinal dis- 
turbances will be deliberately introduced into the chamber by firing an 
external powder pulse to enable a study of their growth and decay under 
different environmental conditions and to shed light upon previously ob- 
served spontaneous waves. 

A perturbation analysis will be carried out on the steady-state model 
to obtain an expression for the acoustic admittance of the burning turbu- 
lent boundary. These theoretical studies will also be supported by addi- 
tional experiments in the slab burner and the wind tunnel combustion simu- 
lator . 
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Appendix A 


DIFFUSION-LIMITED THEORY OF STEADY-STATE HYBRID COMBUSTION 


In the steady state, all of the heat transferred to the solid is used 
in vaporizing the solid material. Under these conditions the heat trans- 
ferred to the surface is balanced by the mass flow rate of fuel times the 
total heat of gasification: 


Q = m h = p r h 
w f v f v 


(A-l) 


The heat of gasification h is composed of the latent heat of vaporization 

v 

plus the heat required to raise the temperature of the solid from ambient 
conditions deep in the grain up to the surface temperature. 


The total heat transfer is generally composed of convective and radi- 
ative parts. The convective heat transfer can be directly characterized 
by the semiempirical theories which are available, whereas a description 
of the radiative transfer relies currently upon an empirical radiation 
coefficient, as will be discussed below. 


Convective Heat Transfer 

The convective heat transfer at the surface is governed by the con- 
duction equation: 

( k 9h\ 

— i7) -W <A - 2) 

P /w 

where the subscript b refers to conditions at the flame. This equation 

serves to define a Stanton number C which is based upon the sensible 

H 

enthalpy difference between the flame and the wall [Ah = (c T) - (c T) ] . 

p b p w 
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Equations A-l and A-2 can be combined to give an equation for the regres- 
sion rate: 


P 


r = 
f 


Ah 

C p u — 
H K b b h 

v 


(A-3) 


It is convenient to evaluate the Stanton number in terms of the skin 
friction coefficient for which experimental data and we 11 -documented em- 
pirical expressions are available. For a nonunity Prandtl number but a 
unity Lewis number, the heat transfer in the boundary layer is related to 
the shear stress by 


Q 


ah/^y au/ay 


T -0.67 

Pr 


(A-4) 


Since Q/Q = t/t from Eq . 4, for a constant Prandtl number Eq . A-4 may 
w w 

be integrated between the surface and the combustion zone to yield 


Q T 

w w -0.67 
— = — Pr 
Ah u 

b 


(A-5) 


Dividing by p u and using the usual definition of the skin friction co- 
b b 2 

efficient, c /2 = t /p u C is related to c /2 by 
f w r e e H f 


c „ n u 

f -0.67 M e e 

C = — Pr 

H 2 2 

p u 
H b b 


(A-6) 


If, in addition, a mass transfer number B is defined by 


u 

e Ah 

B “ u h 
b v 


(A-7) 


the regression rate Eq . A-3 becomes for a unity Lewis number: 


* The present treatment may be easily extended to the case of Le 4 1 
Appropriate corrections are discussed by Lees. 14 
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r = 


B 


(A-8) 


P 


f 


c 
J 
"e e 2 


f -0.67 
p u — Pr 


The mass transfer number B may be regarded as a thermodynamic con- 
stant that characterizes the propellants, since the value of B is deter- 
mined by physical properties of the gas-phase and solid-phase components. 
Also, since B 7 = Pr ' B is the similarity parameter of a boundary layer 
with mass injection, the velocity profiles will have a similar shape every- 
where in the boundary layer if B / is constant. 1 In a turbulent boundary 
layer, Pr is very near to unity so that B 7 and B are not just proportional 
but equal for practical purposes. Thus, B really plays a dual role as a 
thermochemical parameter and as a boundary layer flow parameter. 


To put Eq. A-l in a useful form, it is necessary to evaluate the skin 
friction coefficient c^ for a turbulent boundary layer with mass injection 
and combustion. The skin friction coefficient of an incompressible tur- 
bulent boundary layer is given by 21 


c 

f 

— - = 0.03 Re " ‘ (A-9) 

2 x 

It has been shown 1 > 22 that the effect of mass injection upon the friction 
coefficient can be expressed by 


c 


f 


c 

f 


o 


Jtn(l + B 7 ) 

0.8 

r , „ , 2 

13B 4B 

1 + 10 + 11 

B 7 


B 7 2 



(1 + B 7 ) (1 + 


0.2 


(A-10) 


In the range 5 < B 7 < 100, which corresponds to typical hybrid operation, 
Eq . 10 is adequately represented by the simpler formula: 


f ,-0.77 

= 1.2B 


(A-ll) 
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The presence of combustion manifests itself upon the skin friction coef- 
ficient by producing a variable property boundary layer. Variable property 
effects have been analyzed by many investigators, including Mager, 23 
Burggraf, 24 and Crocco, 26 Marxman 6 has shown that the resulting skin 

friction correction appears as a multiplying density ratio of the form 

- 0.6 

(p/p ) , where p is a suitably defined reference density. For most 

hybrid propellants of interest p/p can be taken as unity to a good ap- 
proximation , s 

The final expression for the regression rate with convective heat 
transfer alone is 

- 0.6 ,0.23 -0.2 

p r = 0.036G(p/p ) B Re (A-12) 

i e x 

A notable fact revealed by Eq . A-12 is the strong coupling between the 
convective heat transfer and the mass injection rate which manifests it- 
self in the small power on B / . This implies that large changes in either 

the heat of reaction (i .e . , Ah) or the heat of gasification h will only 

v 

induce small changes in the regression rate. 


Radiative Heat Transfer 


The heat transfer by radiation is formulated by analogy with an 
equilibrium flow, even though nonequilibrium conditions certainly exist 
in the hybrid boundary layer. This is done because little is known about 
nonequilibrium radiation. In this case 


Q = a e e ( T 
r w g \ r 



(A- 13) 


where the parameter e represents the emissivity of the radiating field. 

g 4 4 

In all practical cases T « T and can therefore be neglected. A new 

w r 

formulation for e is given in Appendix B. 
g 
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It has been shown 6 that the convective and radiative heat transfer 
cannot be directly added because of the coupling that exists between the 
convective heat transfer and the mass injection rate. The energy of radia- 
tion tends to increase the rate of fuel vaporization, and this increase, 
in turn, tends to decrease the rate of convective heat transfer to the 
wall. The resulting trade-off can be expressed by 6 



(A-14) 


where Q is the convective heat transfer that would be present in the ab- 
c 

sence of radiation. For small values of Q /Q (such as would occur in a 

r c 

typical hybrid with low metal loading in the grain) , the trade-off between 
the terms is nearly exact, and r can be calculated with little error by 
using Q alone . 


Treatment of Particle Loading in the Grain 

The preceding section considered the effect of radiation on the hybrid 
burning rate. Many systems that contain an appreciable metal loading in 
the fuel exhibit a significant radiative heat transfer to the grain be- 
cause of the presence of solid combustion products. Systems also exist 
that have a high particle loading in the surface but exhibit little radi- 
ation because the products of combustion are gaseous. Therefore, it is 
necessary to determine the effect of particle loading on the regression 
rate of both radiative and nonradiative systems. 

It can be shown that the presence of particles in the flow has a 

negligible direct effect on the shear stress distribution or velocity 

profile shape in the boundary layer . 26 From this, one can infer that the 

variation of c /c with surface injection should also be unaffected to 
*o 

first order by the addition of particles. However, it is physically 
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apparent that the blocking effect, which reduces c , should depend upon 
only the gas mass addition rate at the surface so long as the particles 
occupy negligible volume compared with that occupied by the gas, as they 
do in any practical application. Based on this assumption, the regression 
rate equation should be written in terms of the bulk density of the vapor- 
izing component. This can be accomplished by defining an effective heat 
of gasification by the equation: 


p h = p h = p 

V V f V V 

eff 



= P 


v 



K 

1-K 


c AT 
m 


(A-15) 


This equation describes the fact that the one component is both heated to 

the surface temperature and vaporized, while the other component is only 

heated. The increase in the effective value of h with increasing metal 

v 

loading tends to counteract the increase in the surface heat transfer, 
which is caused by the surface mass addition. These results indicate that 
the regression rate of a grain can be significantly increased by the mere 
addition of particles whether or not these particles provide any signifi- 
cant radiative heat transfer to the grain. The important grain density 
from the standpoint of regression rate is the bulk density of the binder. 


Evaluation of the Mass Transfer Number 

The mass transfer number B which is a governing parameter in the con- 
vective heat transfer equation, can be evaluated according to the scheme 
of Marxman. 6 This requires the calculation of both the sensible enthalpy 
difference between the flame zone and the regressing surface and the ef- 
fective heat of gasification. The latter quantity is given by Eq , A-15. 

It can be shown 21 that the sensible enthalpy at the flame zone is equal to 
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the heat of reaction per unit mass of propellant (fuel and oxidizer) at 
the local mixture ratio which may be taken as three-fourths of the stoi- 
chiometric value . 


The heat of reaction per unit mass of reactant is given by 


Q 


c 


Ah 


Ev Ah° - Ev Ah° 
r r p p 

” Ev M 
r r 


(A-16) 


A sample calculation for the Plexiglas*-oxygen hybrid system would proceed 
as follows: The stoichiometric mixture ratio for this system is 1.92. 

Assuming that the local 0/F ratio at the flame is 1.4, the balanced reaction 
is 


C H O + 4.37 0 - 2.26 CO + 2.74 CO + H 0 

5 8 2 2 2 2 

Then, using standard heats of formation, the heat of reaction per unit 
mass is 


Q 

c 


1 (-92 .8) 


2 .26 (-26 . 4) - 2 , 74 (-94 , 1) - 4(-57.8) 
100 + 4.37 (32) 


3 

10 


= 1900 cal/g = 3430 Btu/lb 


(A-17) 


Assuming the wall temperature is 600°K and the ambient temperature deep 
in the grain is 300°K, and using 0.4 cal/g -°K as the specific heat of 
the species at the wall , 


Ah = 1900 - 0.4(300) = 1780 cal/g (A-18) 

Using a value of h of 350 cal/g results in 



v 


* Trademark , Rohm & Haas , Co . 
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The mass transfer number B is now related to Ah/h by the curve in 

v 

Fig. A-l (which is reproduced from Ref. 27). For the example above, 
the result is B = 9.3. 



FIG. A-1 THEORETICAL HYBRID MASS TRANSFER NUMBER 
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Appendix B 


GAS- PHASE EMISSIVITY 

A calculation of gas-phase emissivity can be carried out by begin- 
ning with a fundamental description of the radiative process. Consider 
the radiation delivered to a cylindrical wall from an annular zone of 
high-temperature particles. In a hybrid motor with a cool core the radia- 
tive heat transfer comes from this type of annular region. The radiation 
delivered to an incremental area dA lying in the x-z plane (see Fig. B-l) 
from volume dV is 

I dV 

d $ = T sin 3 cos cp dA (B-l) 

4nr 

where I is the intensity generated per unit volume at the location of dV 
and t is the transmissivity of the intervening medium between dV and dA. 



(a) GENERAL SPHERICAL COORDINATES (b) ORIENTATION OF RADIATING ANNULUS 

TA- 6638 - 5 


FIG. B-1 COORDINATE SYSTEM FOR CALCULATION OF EMISSIVITY 
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Assuming the emissivity of individual particles to be unity, 


2 4 

I = 4nr ncrT (B-2) 

p r 

The transmissivity is obtained by assuming the absorption process to be 
of first order; i .e . , 

dl 

— = yi (B-3) 

dr 


where y is a constant. The solution is 


I -vr 

T = — = e 

I 

0 


(B-4) 


A photon of light will be intercepted by any particle whose center lies 


within the distance r from the path of the photon. If the particle emis- 

P 

sivity (absorptivity) is unity, the absorptivity coefficient y is equal 
to the average number of particles encountered per unit path length: 


2 

y = rtr n (B-5) 

P 


Using 


dV = r sin (3 dr dcp d[3 


(B-6) 


as the differential volume in spherical coordinates, the above equations 
can be combined to give 


2 

d $ 


2 

cnr 


P 


4 

T 

r 


2 

-rmr r 


e 


P 


2 

sin 


(3 cos cp dr dcp dp dA 


(B-7) 


The gas-phase emissivity is then obtained by integrating over the volume 

4 

after dividing by crT : 

r 


e 

g 



JJI 

P cp r 


2 

-rmr r 


e 


P 


2 

sin 


P cos cp dr dcp dp 


(B-8) 
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The integration limits for the integrals of Eq . B-8 can be established 

as follows: For cp > cp = sin 1 (R /R) , the r integration can be carried 

e 0 

out directly between r = 0 and the value of r established by the equation 


2R cos cp = r sin [3 


(B-9) 


For cp < cp the r integration must be split into two parts to account for 
e 

the fact that no radiation is generated in the cool core. The two regimes 


of r integration are given by 


0 < r < s = 
1 


2 2 2 — 
R cos cp - (R^ - R sin cp) s 


sin (3 


2R cos cp 2R cos cp 

l _ s < r < *■ = s 

sin p 1 sin 3 2 


(B-10) 


The integration over 3 proceeds from a lower limit of zero to an upper 
limit which is dependent upon the z-coordinate of the end of the annulus 
and upon the value of the angle cp. At the end of the annulus 


3 = 


-l 

tan 


2R cos cp 

-*■ = 3 

z e 

e 


(B-ll) 


The final integral must then be multiplied by four to account for all 
four quadrants. Using the results above, 


. cp p 
4 p e p e 


e = — 
g tt 


I J 

0 0 


2 2 2 
-nnr s -nnr (s - s ) -nnr s 

pi P 2 1 p 2 

1 - e + e - e 


4 aTT/ 2 , e 

TT " cp *1 0 L 
e 


1 - e 


2 

-rmr s 

p 2 


cos cp sin 3 ^3 dcp 
(B-12) 


cos cp sin 3 d3 dcp 


Because of the nature of the exponential functions in Eq . B-12 the inte- 
grals must be evaluated numerically. 

The particle number density n that appears in this analysis can be 
calculated from 37 
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n (x) 


3\ K 
4n 1 - K 


1 


(B-13) 


P r g 
P P 


V pX 

t— p r P dx 
m 0 v 


for a grain containing nonvolatile material, and 


3K o 

1 1 M e n x 

n (x) = — I p r P 

4tt 3 m ^0 K v 


dx 


p r g 
P P 


(B-14) 


for a completely vaporizing grain. The density p is the bulk density of 

r v 

the vaporizing component of the grain. 


An expression similar to Eq . B-12 can be developed to account for 
gas-phase radiation in systems that do not produce particulate combustion 
products. In this case the absorptivity coefficient y is proportional to 
the pressure 


Y = aP 

where the proportionality constant a must be obtained experimentally. 
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Appendix C 


DIFFUSION-LIMITED REGRESSION RATE CALCULATIONS 


A computer program based on the work in Appendices A and B was written 
in BASIC language for solution at a time-sharing computer terminal. A 
complete description of the program, including a sample output format, is 
given at the end of this discussion. The regression behavior of circular 
cylinders is described, such as those studied experimentally during the 
current program. 

Theoretical calculations of the mass transfer number for the poly- 
urethane (PU) and polybutadiene-acrylic-nitrile (PBAN) binder systems have 
been carried out as a function of metal loading, with oxygen considered 
as the oxidizer. The results are shown in Fig. C-l . 

The gas flow rate m , which enters both the convective heat transfer 

g 

expression and the expression for the particle number density, must account 



0 0.1 0.2 0.3 0.4 0.5 


ALUMINUM MASS FRACTION 

TA-6638-2 

FIG. C-1 VARIATION OF MASS TRANSFER NUMBER 
WITH ALUMINUM LOADING FOR PU 
AND PBAN BINDERS 
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for the oxidizer used in the production of particles. It is composed of 
the head-end oxidizer flow rate plus the gas evolved from the surface 
minus the oxidizer used to form the combustion particles. The resulting 
expressions are 27 

m = m + ( 1 - — ) P p rp dx (C-l) 

g o \ 1-K/ J 0 H v 

if particles in the grain produce particulate combustion products, 


. . 1 f.X 

m = m + p rp dx (C-2) 

g o 1-K 0 v 


if particles in the grain produce only gaseous combustion products, and 


m = m + (1 - K - K C ) f X p rp dx (C-3) 

g o 1 1 1 J 0 v 


if a completely vaporizing grain produces particulate combustion products. 
The mass of oxidizer consumed in producing the product particles is ac- 
counted for by £ and Q . 

The results of regression rate calculations for the hybrid systems 
utilized during the current program are shown and discussed in Section IV, 
Experimental Studies. It should be emphasized that the analysis embodied 
by the computer code described below is restricted to the transport-limited 
domain of hybrid combustion wherein chemical kinetics do not exert any 
controlling influence. 


The inputs required for operation of the program are shown in Fig. 
C-2. The equations which are solved can be concisely stated as follows 
(Cases 1, 2, and 3 referred to below are described in Fig. C-2): 


p r = p r (1 - K) 
v f 


eff 




(C-4) 


78 



WY-FXP 


12:35 


SF MON 09/16/68 


1 00 

101 
102 
1 03 
1 0-a 

105 

106 

107 

108 
109 
I 10 
I 1 1 
1 12 
1 1 3 
1 1 4 
1 1 5 
1 1 6 
1 1 7 
I 1 8 
1 1 9 
19Q 
1 9 1 
1 22 
1 23 
1 94 

1 25 
1 26 
1 97 
129 
1 30 
1 3 1 
1 32 
1 3 3 
1 3 4 
1 35 
1 36 
1 37 
1 36 
1 39 
1 40 

1 4 1 
1 4P 

1 43 
1 4 4 
1 45 
1 46 
1 47 
1 48 
149 
1 50 


THF PP0GRAM FNTITLFD "HY" IS BASFD UPON HFAT-’t RANSFER 
LIMITFD HYRRID THF0RY. IT CALCULATES THF RFGPFSSI0N BE- 
HAVI0R 0F A CIRCULAR TUBF FUFL CRAIN. THF INPUTS WHICH ARF 
NFCFSSARY F0R OPERATION OF THE PROGRAM ARF DFFINFD AS 
FOLLOWS ACCORDING TO THE LINF NUMBERS GF "HY": 

50 DO--INITIAL INNFR DIAM. (IN.) 

51 L --LFNGTH (IN.) 

52 D2--3UTFR DIAM. (IN.) 

53 K --MASS FRACTION OF NON-VGLAT ILE SURFACE MAT 'L* 

84 K 1 - -MASS FRACTION 0F VOLATILE SURFACE MAT *L. 

WHICH FORES PARTICULATE CHMHUSTI0N PRODUCTS 

55 K2--T0TAL DFNSITY OF FUFL CRAIN (LB/IN3) 

56 K3--FFFFCT IVE HFAT 0F GASIFICATION (BTU/LB) 

57 MO--0X1DIZER FLOW KATF (LB/SFC) 

(IF MO VARIFS, IT MAY FITHFR Lit READ AS 
DATA 0P DFSCRIBED BY AN FOUATION.) 

58 PO--CHAMRFR PRFSSURF (PSIA) ( SFF THF AB0VF . ) 

59 J --CASF NUMBFR 

J=1 IF PART ICLFS IN THF GRAIN PRODUCE 
PART ICI'LATF COMBUSTION PRODUCTS 
J=2 IF PAPTICLFS IN THF GRAIN PK0DUCF 
GASFSUS COMBUSTION PRODUCTS 
J=3 IF A CCMPLF T FLY V0LATILF GRAIN PR0DUCFS 
PART ICULATF COMBUSTION PRODUCTS 

60 UO--0XIDIZFR TFMP . AT INLFT (OP) 

61 I !5- -COMRI 1 ST I ON TFMP. AT AN 0/F RATIO WHICH IS 
3/4 0F THF STOICHIOMETRIC VALUE AT 5 AIM. 

PRFSSURF (OR) 

64 C7--0AS-PPASF RADIATION CONST. 

65 K4- -MASS 0F OXIDIZER CONSUMF D rFR UNIT MASS OF 

NCN-V0LAT ILF SURFACE Mi A T ' L • IN PRODUCING 
PART ICULATF C0M;PUSTI0N PRODUCTS 

66 K5--MASS 0F OXIDIZER CONSUMF'F' PFR UNIT MASS OF 

V0LAT I LF SUPFACF MAT'L. IN PRODUCING 
PART ICULATF COMBUSTION PRODUCTS 

67 CB--MASS OF PARTICIILATF COMBUSTION PRODUCTS 

F0PMFD PFR UNIT MASS OF NON-VOLATILE 
SURFACE MAT'L. 

A9 FP--FMISSIVITY (ABSORPTIVITY) OF FUFL SURFACE 
70 P --MASS TRANSFFP NUMBFR 

73 K6--DFNSITY OF PART ICULATF COMBUSTION PRODUCTS (LB/IN3) 
77 C2--B0UNDARY LAYER DENSITY CORRECTION FACT0K 
76 S --CXIDIZFR VISCOSITY AT TEMP. 110 (LD/IN-SFC) 

82 T1--TIMF INCRFMFNT F0R CALCULATIONS (SFC) 

63 T2--QUIT TIMF (SFC) 

84 T3--TIMF INCREMENT F0R PRINTOUT (SFC) 

85 X1--DISTANCF INCRFMFNT FOP CALCULATIONS (X/L) 

103 PUT IN THF NAM.FS 0F THE FUFL AND OXIDIZER BY 

TYPING F.G.. 

103 PRINT "PMM"#"OXYGFN" 


FIG. C-2 SUMMARY OF INPUT REQUIRED FOR HYBRID COMPUTER CODE 
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rj l>^ 


Q 


eff 


= 0.036 



(C-5) 


m = m + ( 1 - — \ r p i 
g ° \ l-K J J H v 


r rr D dx (Case 1) 


■ . l r> . 

m = m + p r tt D dx 

g o l-K J v 
o 


(Case 2) 


(C- 6 ) 


m + m 


ti + (l-K - K H f n r n D dx (Case 3) 
o 1 1 1 J v 


TT *2 

A = - (D - 26 ) 

P 4 


(C-7) 


* 

26 = 


0.8 

0.21 D (x/5D) x £ 5D 

0.21 D x > 5D 


(C- 8 ) 


Q = a e e T 
r w g r 


(C-9) 


cp P 
e e 


r p I ~ YS 1 “ y(S 2~ S 1 ) _YS 2| 2 

I 11 - e +e -e f cos 9 sin (3 dp dcp 


o o 


tt/ 2 M e f -ys 

J I 

cp o 
e 


2 \ 2 
e I cos cp sin p dp dcp 


(C-10) 


Y = 


tt nr (Cases 1 and 3) 

P 

at p (Case 2) 


(C-ll) 
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(C-12) 



(C-13 ) 

(C-14) 

(C-15) 

(C-16) 


4. Calculate Q 

r 

5. Calculate r 

The details of this computational procedure are shown in the printout of 
the computer code given in Fig. C-3 . 
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PY l r> : | 7 OF MIN 09/1 6/6F 

A? I.FT F3 = 0 
S('LFTnn = l 

b 1 l.r 1 L = 1 o 

SO|_F th? = o 

b3 LFT K = n 

S/: LFT K1=0 

bb I.FT KP=3 • FO 1 06F-P 

Sft I.F T X 3 = ft 3 n 

S 7 I.F 1 PO = .OF. ft 

sf i.ft pn=pnn 

5 9 ITT J = P 

60 I.FT 1 ' 0 = 1 6 0 

6 1 I FT Lb = 3 9.3 6*1 .6 

6?1.FTLP = LS + ( 1 +3 .33F - 6* (PO-7S ) > 

63 LFT F 6=3 
6 4 LFT r 7-0.00 3 
6 b LFT K4 = 0 . RFK 
66 LFT K S =0 

6 7 LFT 0F = 1.f>f 

6F LFT Fl=3.307F-1b 

69 LFT FP=.9 

70 LFT R = 9.ft 

7 1 LIP 0(1 Of. ) 

7P IF IN N(iori) 

73 ITT FR = 0. |437 

7 A LFT C6=(3*9.5AF 4)/(4*Kb*P6) 

7 S LFT J 1 = 1 

76 LFT 01=0.036 

77 LFT CP=1 

7 F LFT S=1.0P F-6 
7 9 LFT 03=(SFL 7 to.o 
FO LFT 0 S= 3 . 1 A) S 9 
F I LFT OF=Cb/Fi 
FO LFT T I =0.05 


F 3 

LFT TP = 11) 



FA 

I.FT T 3 = • b 



Pb 

LFT X |=0.0o 



F 7 

LFT L=0 



8F 

LFT III =0+1 '0/3 



90 

Rp I NT 



96 

LF T 14 = INT (VI /X 1 +1 /? ) 
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POINT 
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1P6 rFINT "f' FTOL POT T . rOT • " » "4F'I A">"idF To- 1" .* "L*F 6 I'A" 

]PH PR I NT K,K|,K4 # KS ,06 
I 99 PRINT 

130 POINT "I 1 FL-T”i , 'nr L-X/L"»"P»IN'l T INF" , " 3L 1 1 1 INF " » "OASF" 

139 PRINT T 1 , X 1 » T3 » TP » J 
134 POINT 

1.35 PRINT 

140 LFT 06 = 01 *r,9*C3*<R> «0.P,3 
145 LFT T4=0 

147 LFT T5=0 
150 LFT T =0 

155 LFT X = X1 

156 LFT 1=1 

1 57 LFT 11=14 
160 LFT N = NC 

165 IF T>0 THFN 175 

170 L.FT nci) = no 

175 LFT R.3=(X*L)/C5*D< I ) ) 

1 HO IF P 3 < = I THFN 1 90 

1 87 IF P3>5 THFN POP 

IBP IF R.3 > 1 THFN P10 

190 LFT F = O ( I ) * ( 1 - 0 • P 1 *(R3)tO.B) 

195 LFT IJ=HO 

1 96 GO T0 PPO 

POO LFT F=Of I )*C] -0.P1 ) 

P05 LFT 1 1 = 1 !P 

PO 6 GO TO PPO 

RIO LFT F=0 < I )*< 1 -C .9 1 ) 

PI 5 LFT U = HO + ( (0?-HC)/4) + (i,'.3-l ) 

PPO LFT 4=04*(F)tP 

PP 5 LFT 01=06*(X)» <-O.P)*CN/A)tO.H 
P.30 LFT 0.3 = 05*01 *DO*XI *L 
P.35 LFT PI =01 /(KP*f 1 -K) ) 

P 40 IF .1=1 THFN P55 

P45 IF J=P THFN 345 

P 50 IF J = 3 THFN 300 

9 5 5 IF I>1 THFN P PO 

P60 LFT N = NO + ( 1 - K *K 4 / ( 1 -K) ) +°3 

P61 LFT N1 =160+03*1/(1 -K> 

P 6 5 LFT 05 = 03 

P70 LFT N C I >=06*CF +K/C 1 -K)*PC/I'*1 /N *03*(P3>* C-Q.h) 

P71 GO T0 366 

P80 LFT N = N + ( 1 - K * K 4 / ( I - K ) ) * 0 4 
PHI LFT M =M +04*1 /d -K) 

PB5 IF R3 > 1 THFN P95 

P 90 LFT N C I)=06*08*K/C 1 -K ) *rO/l ! *1 /N *06* !KI) t ( -o . H) 

P 9 1 GO T0 3 66 

295 LFT N( I )=06*08+K/( 1 -K ) *P0/l i* 1 /N +0 5 

P 96 GO T 0 366 

300 IF I>1 THFN 3P5 

305 LFT N=NO+( 1 -K1 -K1 *K5)*03 

306 LFT 16 1=^.0+0 3 
310 l.FT n 5 = 3 3 

315 LFT N { I )=C6*K1 *l'0/li*l /X*».3*(F\3)t (-0.B) 

316 GO TO 366 

3P5 LFT N=N+< 1 -K1 -K1 *K5)*04 
396 LFT M =<61+0 4 
330 IF P3>1 THFN 340 

335 LFT N( I >=06*KI *H0/H*1 /N*'15* (P.3 ) t C -0 . H ) 

336 GO TO 366 

340 l.FT M I ) = 0 6 *K 1 *P0 /1 1 * 1 /N *15 


FIG. C-3 Continued 
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3 4 1 GG 1 r 366 

345 IP I>| THFN 360 

360 LFT ,v = + C ] / ( ] - K) ) > 3 

36! LFT r-l=F 

355 LF'l N ( I ) = 0 7 + F ( 1 

366 on TO 366 

360 LFT X=i' + ( 1 /(] -K) 

36! LFT N1=X 

365 LFT Ntn = 07FP0 

366 IF J1=0 TFFN 371 

36 3 ITT ft> = pi 
370 OB TO 390 

37 1 IF J = P THFN 354 

3 7 •3 IF i\l(I)>0 THFN 375 

373 LFT PP = «I 

374 GO TO 390 

375 IF I <14 THFN 3 HP 

376 IF I >6+1 4 THFN 3 HP 
37 F IF 1=1! THFN 3 HP 

9 HO 0 0 TO 3! 7 

3H° IF 1 = 15 IPIN 590 

303 PC 10 3 ; 7 

3F/i IF I >”+ f4 -|HF\ 36 7 

3HS GO TO 375 

3 F 7 l.FT 1P = M *F9/K3*F3 + C" 1 ) t 4 

•7F H LFT FO = Of + ( T’/ 1 1 +F ,\P C - )P/ ; 1 ) )/(KP* ( 1 -K) ) 

3 90 LFT O4=n5+K0*( 1 -h )+i'P+n( j )*xi *L 
391 LFT 05 = 115+04 

395 IF 1=11 TFFiN 495 

400 I.FJ n ( I ) = n< I ) +P+FP + T 1 

401 IF n(I)>OP THFN 795 
405 IF X»0 .995 THFN 4?5 
410 LFT 1=1+1 

415 LFT X=X+X1 

4P0 GO TO 165 

40 5 IF T = 0 TFFN 450 

430 l.FT P3= CX 1 +H0 3/P-H0 

435 LFT '.■'=V+V3*T1 

440 LFT i v P = M 

440 IF T = T 4 THFN 455 

444 LFT T = T +’M 

446 IF T > T 0 1FFN HOO 

4 4 F GO TO 155 
460 LFT Ki> = n 

450 IF J1=0 THFN 455 

453 LFT J1=P 

454 GO TO 1 55 
456 FLINT 

460 HP I. NT "oor:. I 'F L- F" , " l- T . LOSS" »" GAS FLOV", "TOTAL FLHV " , "G-0" 

46! l.FT ri6 = ('.■■/ (04*K°+L ) +CPO ) tP ) »(). 5 

46 5 r F INI (06 -00 ) /O ,'-'*4 S3 . 6*N»f- 1 , 6 0 / C 0 4 + C I >6 ) » 2 ) 

470 HP I NT 

471 P P I NT 

475 LFT T = T +T 1 

476 l.FT T4 = T4+T3 

4 1 -’0 IF T >T0 THFV HOP 
490 GO TO 155 
4 95 IF' T= 1 4 TFFIM 6ii6 
496 I FI 0(1 ) - O ( J ) + r Ji' pt'll 
600 I H 11 = 11+14 

505 PO TF 40 6 

506 IF .11 =? T H F \ 5 1,0 


FIG. C-3 Continued 
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50K Gn TP 405 

510 IF 11 = 14 THFN 520 

515 00 T0 535 

520 PRINT "T IMF" 

525 PRINT T ' 

526 PPINT 

530 PRINT "niSTANCF'V'RFG. HATF ,, «"DFL-K"«"ARFA ,, *"FMIS*I VJTY" 

531 PRINT 

535 PPINT X,R2,(r'(I)-nO>/2.C4*<nCI>>»S?#F3 

536 PRINT PI 04*(F>»« 

540 LFT r>< I > = 0< I >+2*RP*TI 

54) ip r>(I>»n? THFN 795 
545 LFT 11=11+14 
550 GO TO 405 
590 LFT TS=T5 ♦ (?*T3 ) 

600 LFT R=n<I)/P 

605LFTR0=R-0.P5 

610 LFT Y1 =3. 1 41 59/30 

615 LFT Y2 = ATN< 1 /< (R*R 5/(K0*R0>- 1 ) t (0.5 > ) 

60S LFT F3=0 
63 1 LFT J5 = 1 

633 LFT K9= 1 

634 LFT Y= Y 1 / 1 00 

635 LFT 2=3.14159/0 
640 LFT J4= 1 

645 LFT Y3=?*P*r0SCY) 

646 IF Y> YP THFN 650 

648 LFT 53= < Y3/2-(R0*R0-R*R*(SIN( Y> ) *2) » (0.5 > )/SIN( 2) 

650 LFT X3=J4*XI*L 
655 LFT 54=Y3/SIN(2) 

660 IF Y<Y? THFN 666 
66? LFT K9=0 
664 LFT S3=0 

666 IF J5= 1 THFN 670 

667 LFT S5=K9+FXP(-N(I-J4)*S3>-K9*FXP(-N< I - J4 ) *( S4- S3 ) )+EXP< -N( I- J4>*S4) 

66o GO TO 675 

6 70 LFT 5 5 = K 9*F XP ( - N ( I +J4 > *53 ) - K9*F XH ( -N ( I ♦ J4 > * ( S4-S3 > > +FXF ( -N ( I+J4>*S4) 
675 LFT 56= (1 -S5)*(.5IN< Z) )»2*CQ5( Y> 

680 LFT Zl=(Y3*Xl*L)/((X3>t2+(Y3>»2> 

685 LFT 57=56*21 *Y1 *?/3. 1 41 59 
690 LFT F3=F3+57/l .06988 
695 LFl 2=2-21 

700 LFT J4 = J;+1 

701 IF J5=l THFN 705 
70? IF J4>I THFN 7?0 

703 IF J4 > 1 0 THFN 720 

704 GO jo 645 

70 5 IF J4>l/M-I THFN 720 

706 IF J 4 > 1 0 THFN 720 

715 GO TO 645 

7?o LFT Y = Y + Y 1 

7?5 IF Y>3. 14159/2 THEN 745 

740 GO TO 635 

745 IF J5=? THFN 750 

746 LFT .15 = 2 
748 GO TO 633 
750 GO TO 387 

795 PPINT " 5PFG IF IFF) FXTFPNAL H I h8 F IF « F XGF FDFP" 

801! FNO 


FIG. C-3 Concluded 
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Results of a typical calculation are shown in Fig. C-4 for a 1-in. 
Plexiglas (polymethylmethacrylate) grain. At each position, the regres- 
sion rate that would have been obtained in the absence of radiation is 
printed out below the total regression rate. In this particular system, 
the radiation is so low that it produces little effect on the regression 
rate, partly because of the trade-off between radiative and convective 
heat transfer. Two port areas are also printed out at each position; the 
first is the geometric area and the second is the true flow area that is 
the geometric area minus the area occupied by the boundary layer displace- 
ment thickness . The flow area must be used in the calculation of the mass 
flux which enters the convective heat transfer equation. Results obtained 
for other hybrid propellant combinations are given in the main text of this 
report . 
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Appendix D 


DEVELOPMENT OF A MODEL FOR THE PRESSURE-SENSITIVE DOMAIN 


At high pressures a thin turbulent diffusion flame is established in 
the boundary layer at a position where the local mixture ratio will allow 
combustion. Under these conditions the gas-phase reaction rates are large 
compared to the turbulent transport rates, and the latter are then the 
rate-controlling factors. As the pressure is decreased, the reaction 
rates decrease, the flame zone broadens, and the combustion process ap- 
proaches that of a premixed turbulent flame; i.e., the combustion time 
becomes long compared to the mixing time. 

An important point to consider is whether the gross kinetics effect 
on the regression rate is caused by the behavior of the heterogeneous sur- 
face reactions associated with the decomposition of the solid or by the 
behavior of the gas-phase reactions. When gas-phase reaction rates are 
high, the mainstream oxidizer mass fraction must be low between the flame 
zone and the wall, because high reaction rates imply a thin flame zone and 
consequent effective utilization of the oxidizing species. Thus, in this 
case, any heterogeneous reactions at the surface will be controlled by the 
diffusion of the small proportion of oxidizing species that is available 
to the wall. Since, by Reynolds analogy, the diffusional transport of 
mass momentum and energy is governed by the same equation, no correction 
to the basic hybrid model is required to cover this case. It can be con- 
cluded that the presence of pressure-sensitive heterogeneous reactions 
at the gas-solid interface does not explain the observed decrease in re- 
gression rates at low pressure. 
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The foregoing argument indicates that the gas-phase reaction rates 
at low pressure must be examined in order to explain the regression rate 
behavior. The physical explanation of the general process is that at 
high pressures, the gas-phase reaction rates are fast enough to consume 
all of the available material that is being vaporized at the surface. As 
the pressure is lowered, however, the gas-phase reaction rates decrease, 
and finally a point is reached where the flame cannot efficiently consume 
all of the vaporizing material and the local mixture ratio of the flame 
begins to shift. This shift in mixture ratio leads to a decrease in heat 
release, which in turn decreases the vaporization rate of the surface. 

At pressure well below the threshold level where the chemical reaction 
time is much longer than the gas-phase mixing time, the flame behavior 
approaches that of a premixed flame, 

A detailed solution of the complete turbulent boundary layer equations, 
involving not only chemical kinetics but also the interaction between the 
turbulence mechanism and the combustion process, is obviously beyond the 
current scope of our understanding. These considerations naturally lead 
to an inquiry into whether the simplified turbulent flame theory developed 
in the past for premixed flames 28 can be applied to the hybrid case in 
which there is a transition from a premixed flame to a diffusion flame as 
the pressure increases. An excellent review of the wrinkled-flame-sheet 
turbulent flame theory has been given by Karlovitz . 28 

The theory of Karlovitz is a one-dimensional theory; therefore its 
applicability to a two-dimensional boundary layer problem must be carefully 
considered. The use of a one-dimensional theory as a first approximation 
is supported by schlieren photographs, which consistently show that most 
of the heat release takes place very chose to the surface even at low 
pressures . 1 Furthermore, in the vicinity of the surface, a one-dimensional 
approximation describes the turbulent transport processes very well , 20 
lending credence to the present approach. 
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Application to the problem at hand leads to a description of the re- 
gression rate behavior in terms of the ratio t/T, where T is the charac- 
teristic mixing time corresponding to turbulent diffusion, and t is the 
characteristic reaction time for a premixed zone of dimensions correspond- 
ing to the turbulence scale. 

The mean-square displacement of a particle from its original position 

28 

by turbulent fluctuations is described by the equation 


t 



o 


If the correlation function is assumed to be given by 


-t/T 

R = e 
t 


Eq . D-l can be evaluated to give the rms displacement: 



The average flame speed in the gas phase is given by: 


(D-D 


(D-2) 


(D-3) 


S = 




T -t /T 

— (1 - e 1 ) 

1 


(D-4) 


In the mixing (heterogeneous) limit, t^/T-O and S-*u ' where u 7 is a charac- 
teristic turbulent velocity in the combustion zone. Since the turbulent 
flame speed must be compatible with the regression rate at the wall , the 
regression rate ratio is given by 


r 

r 

CO 




(D-5) 
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The remaining task is to express the ratio t /T in terms of flow parameters 
such as pressure and temperature. 

From a physical point of view, it is convenient to express the time 
ratio in terms of length and velocity ratios as 


t 


1 


T 


U f *2 


(D-6) 


In Eq . 41, j is the distance that the flame front propagates into unburned 
gas at the (kinetic) flame speed u^ while £ is a characteristic scale of 
turbulence. Using the Denison and Baum formulation 29 for the flame speed 
u and the asymptotic behavior of S to describe u 7 , the time ratio becomes 


t 

_1 

T 


-0.2 0.23 

£ GRe B 

C £ 

£_ n/2 l+n/2 


2 p 


U 


exp(-E f /2RU f ) 


(D-7) 


The constant C is required to complete the kinetic formulation of u^ and 
must be determined by comparison with data. 

The length £ , which can be thought of as related to the integral 
scale of turbulence, is probably nearly constant over a wide range of flow 
conditions. On the other hand, the length £^ could, for example, conceiv- 
ably depend upon the mass flux because of changes induced in the relative 
availability of fuel and oxidizer in the vicinity of the flame zone. by 
changes in the mass flux. At present the disposition of this point must 
be left to experiment. 
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Appendix E 


REGRESSION RATE MEASUREMENTS WITH OXYGEN 

The basic components of the 2-g-in .-diameter motor that was used to 
obtain regression rate measurements in the low-pressure regime are shown 
in exploded view in Fig. E-l . A 12-in, -long flow-straightening section 
is used upstream of the inlet to the grain to ensure uniform flow condi- 
tions. Screens and a porous plate are used at the injector end of the 
straightening section to help spread the incoming stream from the injector. 
At the grain end of the straightening section, an inlet for propane and a 
spark plug are provided as the ignition system. 

An existing flow control facility was adapted for the gaseous oxidizers 

used in the small-scale tests. The facility was designed for liquid engine 

tests with cyrogenic oxidizers and therefore was constructed of 304 stain- 

* + 

less steel with Teflon and/or Kel-F 1 seats and seals. The flow system 
schematic is shown in Fig. E-2 , The oxidizer mass flow is kept constant 
throughout a test by passing the flow through a sonic choke that is oper- 
ated above the critical pressure ratio. 

The grain is ignited by preflowing a small amount of oxygen and pro- 
pane into the chamber and igniting it with a spark plug. The ignition sys- 
tem is preset so that it operates near the lower flammability limit for 
the oxygen-propane mixture, thus yielding an oxygen-rich high-temperature 
source. The time for the flame to spread from the head end of the grain 


* Trademark , E . I . duPont de Nemours and Co . 
t Trademark, 3M Company. 
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FIG. E-1 SMALL-SCALE TUBULAR HYBRID MOTOR 


DOME 

LOADER DOME 



FIG. E-2 OXYGEN FLOW-CONTROL SYSTEM 

to the nozzle end was determined to be approximately one second. There- 
fore, the main flow valve was delayed for that length of time. The weight 
of fuel typically consumed during this ignition process was determined to 
be 0.5 g which is about 1 percent of the amount of fuel that would be 
consumed in a 5-sec test at a regression rate of 0.01 in. /sec. 

An automatic sequencer actuates the system. The run duration is set 
and the events proceed as follows; 
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Zero time 


+1.9 sec 

+2.0 sec 

+2.9 sec 

+3 .0 sec 

+8 , 13 , or 
18 sec 


Pressure-recording oscillograph on 
Oxidizer preflow valve open 
Propane valve open and spark source on 
Oxidizer main flow open 

Oxidizer preflow and propane valves closed 

Oxidizer main flow closed; nitrogen purge 
valve open 


Event markers are recorded on the oscillograph along with pressure-time 
data obtained from strain-gage transducers at the entrance and exit of 
the combustion changer. Running times of 5, 10, and 15 sec were chosen 
to obtain the regression rate/time dependence. 


To check the experimental apparatus, weight loss data were obtained 
for the Plexiglas-oxygen propellant combination and compared with the pre- 
dicted numerical results. This system was chosen because it has been ex- 
tensively documented in the literature. 2 ’ 26 


Data were obtained at 5-, 10- , and 15-sec intervals for chamber 
pressures of 45, 115, and 215 psia. The results are shown in Fig. E-3 . 

The regression rate, and therefore the weight loss, lessens as the pres- 
sure is lowered, indicating that the burning process is not entirely heat- 
transfer controlled at low pressures but that chemical kinetics must play 
a role. Previous results have indicated that the heat-transfer controlled 
limit, corresponding to the theoretical curve, is reached in the neighbor- 
hood of 265 psia 26 for this propellant combination. 
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WEIGHT LOSS 



FIG. E-3 WEIGHT LOSS OF A PLEXIGLAS GRAIN AS A FUNCTION 
OF CHAMBER PRESSURE 
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Appendix F 


CHAMBER PRESSURE RESPONSE TO VARIABLE OXIDIZER MASS FLOW 


The equation expressing conservation of mass in the combustion cham- 
ber is 


d(pV) 

— e — = m 
dt o 


+ m. 


m 

n 


(F-l) 


Using the perfect gas law, p = pRT, and the equation for the mass flow 
through the nozzle , 


m 

n 



(F-2) 


Eq . F-l becomes 


V dp pV dT _ . 

RT dt " 2 dt t 

RT 

where the oxidizer and fuel mass flows have been lumped together for 
convenience. 



Two cases can be easily calculated — an isothermal process and an 
isentropic process. In the isothermal case dT/dt = 0 and Eq . F-3 may be 
written in the simple form 


dP 
*1 dt 


K 


- m 


-V 


(F-4) 


where P = p/p , and 
i 
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(F-5a) 


P.V 

l 

1 ~ RT 

i 


K 


P. 

1 

RT 


Y+l 

Y- 1 


2 

, Y +1 , 


i 


(F-5b) 


Since Eq . F-4 is linear, solutions can be superimposed. A typical forcing 
function of the form applied experimentally is 


m^ = m^ + (m 2 - m^) (1 + sin tut) 


(F-6) 


which denotes the superimposition of a step function and a sine wave 
The solution for P is 


P = 


( m 2 - m^ r 


K. 


-t/T sm((ut + cp) 
1 - e + i 

2 2 i 

(1 + T (D ) E 


+ 1 


(F-7) 


-1 


where cp = - tan Tuu is the phase angle, and t = K /K is the characteris- 

1 2 

tic time constant of the chamber. At large time, as t -» co ; the mag- 
nitude of the pressure oscillation divided by the magnitude of the mass 
flow oscillation is 


R = 


K P 

_2 

m -m 
2 1 


(1 + T 2 tv 2 )' : 


(F-8) 


Y 


In the isentropic case, when p = p , the counterpart of Eq . F-3 is 


K. 


Ill 


_1 p" y dP 

Y 


= M - K P 
dt 2 


Y+l 

2 Y 


(F-9) 


This equation is nonlinear and must be solved numerically. 


100 



Appendix G 


THERMOCHEMICAL CALCULATIONS FOR HYBRID PROPELLANTS 

Thermochemical calculations were carried out for all of the PBAN- 
oxygen and PU-oxygen propellant systems studied experimentally at chamber 
pressures of 5 and 15 atm. Two 0/F ratios were considered, stoichiometric 
and three-quarters of stoichiometric, the latter value most closely repre- 
senting hybrid operation. The calculated flame temperatures are needed 
as inputs to the hybrid computer code . The results of the calculations 
are given on the following pages. 
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HYBRID SYSTEM NUMBER 1A 


PRESSURE « 5 ATMOSPHERES 


ELEMENTS 

IN 

REACTANTS, 

ATOMS 


FUEL 

OXIDIZER 

BINDERS 

H 




1.9756 

10.34 

0 . 

0 . 

C 




1.3279 

6.95 

0 . 

0 . 

0 




4.3061 

0.13 

2.00 

0 . 

N 




0.0401 

0.21 

0 . 

0 . 

AL 




C.4670 

0 . 

0 . 

l.CO 

ENTHALPY 

OF 

REACTANTS, 

KCAL 

-16.3358 





PRESSURE, ATM 
TEMPERATURE, DEG K 
HEAT CAPACITY, CAL 
ENTHALPY, KCAL 
ENTRCPY, CAL / DEG 
FROZEN GAMMA 
SHIFTING GAMMA 
MGLS CF GAS 
MOLECULAR WEIGHT 
DENSITY, GM / CC 


CHAMBER 
5.0000 
3606.9656 
/DEG K 233.5696 
-16.3358 
K 240.5331 

1.1948 
1.1187 
3.3373 
28.028C 
C. C0C506155 


HYBRID 

SYSTEM NUMBER 

1A 

PRESSURE =15 

ATMOSPHERES 



ELEMENTS 

IN 

REACTANTS, 

ATOMS 


FUEL 

OXIDIZER 

BINDERS 

H 




1.9756 

10.34 

0 . 

0 . 

C 




1.3279 

6.95 

0 . 

0 . 

C 




4.3061 

0.13 

2.00 

0 . 

N 




C.C4C1 

0.21 

0 . 

0 . 

AL 




0.4670 

0 . 

0 . 

1.00 

ENTHALPY 

OF 

REACTANTS, 

KCAL 

-16.3358 





PRESSURE, ATM 
TEMPERATURE, DEG K 
HEAT CAPACITY, CAL 
ENTHALPY, KCAL 
ENTROPY, CAL / DEG 
FROZEN GAMMA 
SHIFTING GAMMA 
MOLS CF GAS 
MOLECULAR WEIGHT 
DENSITY, GM / CC 


CHAMBER 
15.0CCC 
3802.0615 
/DEG K 201.1827 
-16.3358 
K 233.3197 

1.1896 
1.1232 
3.27C2 
28.5647 
C. 001470C14 


60% PBAN/40% A1 ( 0/F = 2.16) 
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HYBRID SYSTEM NUMBER IB PRESSURE = 5 ATMOSPHERES 


ELEMENTS IN REACTANTS* 

ATOMS 

FUEL 

OXIDIZER 

BINDERS 

H 

2.3832 

10.34 

0 . 

0 . 

C 

1.60X9 

6.95 

0. 

0 . 

0 

3.9050 

0.13 

2.00 

0 . 

N 

C.C484 

0.21 

0. 

0 . 

AL 

C. 5634 

0 . 

0 . 

1.00 

ENTHALPY OF REACTANTS, 

KCAL -19.7067 





CHAMBER 




PRESSURE, ATM 

5.0000 




TEMPERATURE, DEG K 

3667.0096 




HEAT CAPACITY, CAL /DEG 

K 244.5151 




ENTHALPY, KCAL 

-19.7C67 




ENTROPY, CAL / DEG K 

248.9351 




FROZEN GAMMA 

1.1965 




SHIFTING GAMMA 

1.1198 




MOLS OF GAS 

3.5238 




MOLECULAR WEIGHT 

26.3364 




DENSITY, GM / CC 

G. 000471 5 1C 





HYBRID 

SYSTEM NUMBER 

IB 

PRESSURE =15 

ATMOSPHERES 



ELEMENTS 

IN 

REACTANTS, 

ATOMS 


FUEL 

OXIDIZER 

BINDERS 

H 




2.3832 

10.34 

C. 

0. 

C 




1.6019 

6.95 

0. 

0. 

0 




3.9050 

0.13 

2.00 

0. 

N 




C. 0484 

0.21 

0. 

0. 

AL 




0.5634 

0. 

0. 

l.CO 

ENTHALPY 

OF 

REACTANTS, 

KCAL 

-19.7067 





PRESSURE, ATM 


CHAMBER 

15.0CCC 

TEMPERATURE, DEG K 


3867.8880 

HEAT CAPACITY, CAL 

/DEG 

K 207.98C1 

ENTHALPY, KCAL 


-19. 7067 

ENTROPY, CAL / DEG 

K 

241.3212 

FROZEN GAMMA 


1.1912 

SHIFTING GAMMA 


1.1246 

MOLS CF GAS 


3.4505 

MOLECULAR WEIGHT 


26.8536 

DENSITY, GM / CC 


0.001369447 


60% PBAN/40% A1 ( 0/F = 1.62) 



HYBRID SYSTEM NUMBER 1C PRESSURE * 5 ATMOSPHERES 


ELEMENTS 

IN 

REACTANTS, 

ATOMS 


FUEL 

OXIDIZER 

BINDERS 

H 




2.4041 

10.34 

0. 

0. 

C 




1.4159 

6.95 

0. 

0. 

C 




4.4802 

0.13 

2.00 

0. 

N 




C. 0488 

0.21 

0. 

0. 

AL 




C.2150 

0. 

0. 

1.00 

ENTHALPY 

OF 

REACTANTS, 

KCAL 

-19.8795 





PRESSURE, ATM 


CHAMBER 

5.0CCC 

TEMPERATURE, DEG K 


3430.5849 

HEAT CAPACITY, CAL 

/CEG 

K 230.4661 

ENTHALPY, KCAL 


-19*8795 

ENTROPY, CAL / DEG 

K 

255.0476 

FROZEN GAMMA 


1.2081 

SHIFTING GAMMA 


1.1183 

MOLS OF GAS 


3.6522 

MOLECULAR WEIGHT 


26.6035 

DENSITY, GM / CC 


C.00C486299 


HYBRID SYSTEM NUMBER 

1C PRESSURE *15 

ATMBSPHERES 


ELEMENTS IN REACTANTS, 

ATOMS 

FUEL OXIDIZER 

BINDERS 

H 

2.4041 

10.34 C. 

0 . 

C 

1.6159 

6.95 C. 

0 . 

0 

4.4802 

0.13 2.00 

0 . 

N 

C. 0488 

0.21 C. 

0 . 

AL 

C.2150 

0 . 0 . 

1.00 

ENTHALPY OF REACTANTS, 

KCAL -19.8795 




CHAMBER 



PRESSURE, ATM 

15.0CCC 



TEMPERATURE, DEG K 

3602.9991 



HEAT CAPACITY, CAL /CEG K 199.7885 



ENTHALPY, KCAL 

-19.8795 



ENTROPY, CAL / DEG K 

247.1396 



FROZEN GAMMA 

1.2032 



SHIFTING GAMMA 

1.1232 



MOLS CF GAS 

3.5916 



MOLECULAR WEIGHT 

27.0389 



DENSITY, GM / CC 

C. 001412466 




80% PBAN/20% A1 ( 0/F = 2.54) 
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HYBRID SYSTEM NUMBER ID 


PRESSURE * 5 ATMCSPHERES 


ELEMENTS IN REACTANTS, 

ATOMS 

FUEL 

OXIDIZER 

BINDERS 

H 


2.9268 

10.34 

0. 

0. 

C 


1.9672 

6.95 

0. 

0. 

0 


6.0993 

0.13 

M 

• 

O 

o 

0. 

N 


0.0594 

0.21 

0* 

0. 

AL 


C. 2595 

0. 

0* 

1.00 

ENTHALPY OF REACTANTS, 

KCAL *24.2012 






CHAMBER 




PRESSURE, ATM 


5.0CCC 




TEMPERATURE, DEG K 


3430*6033 




HEAT CAPACITY, CAL 

/DEG K 206*6955 




ENTHALPY, KCAL 


-24* 2012 




ENTROPY, CAL / DEG 

K 

267*0609 




FROZEN GAMMA 


1*2131 




SHIFTING GAMMA 


1*1236 




MOLS OF GAS 


3*9283 




MOLECULAR WEIGHT 


24*6515 




DENSITY, GM / CC 


C. 000452116 





HYBRID 

SYSTEM NUMBER 

ID 

PRESSURE >15 

ATMOSPHERES 



ELEMENTS 

IN 

REACTANTS, 

ATOMS 


FUEL 

OXIDIZER 

BINDERS 

H 




2.9268 

10.34 

0. 

0 . 

C 




1.9672 

6.95 

0. 

0 . 

0 




4.0993 

0.13 

2.00 

0 . 

N 




C.0594 

0.21 

0. 

0 . 

AL 




0.2595 

0 . 

0. 

1.00 

ENTHALPY 

OF 

REACTANTS, 

KCAL 

-24.2012 





PRESSURE, ATM 


CHAMBER 

15.0000 

TEMPERATURE, DEG K 


3596.2933 

HEAT CAPACITY, CAL 

/DEG 

K 175.2233 

ENTHALPY, KCAL 


-24.2C12 

ENTROPY, CAL / DEG 

K 

258.5418 

FROZEN GAMMA 


1. 2085 

SHIFTING GAMMA 


1.1293 

MOLS OF GAS 


3*8664 

MOLECULAR WEIGHT 


25.0932 

DENSITY, GM / CC 


C. 001314573 


80% PBAN/20% A1 (O/F = 1.91) 



HYeRIC 

SYSTEM NUMBER 

IE 

PRESSURE * 5 

ATMOSPHERES 


ELEMENTS 

IN 

REACTANTS, 

ATOMS 


FUEL 

OXIDIZER 

H 




2.5923 

10.34 

0. 

C 




1.7424 

6.95 

0. 

0 




4.7326 

0.13 

2.00 

N 




C.C526 

0.21 

0 . 

ENTHALPY 

OF 

REACTANTS, 

KCAL 

-21.4353 




PRESSURE, ATM 
TEMPERATURE, DEG K 
HEAT CAPACITY, CAL 
ENTHALPY, KCAL 
ENTROPY, CAL / DEG 
FROZEN GAMMA 
SHIFTING GAMMA 
MOLS OF GAS 
MOLECULAR WEIGHT 
DENSITY, GM / CC 


CHAMBER 
5.0CCC 
3286* 74B1 
/DEG K 224.4207 
-21.4153 
K 262.391 2 

1.2174 
1.1172 
3.8294 
26.1136 
C.0004841C1 


HYBRID 

SYSTEM NUMBER 

IE 

PRESSURE *15 

ATMOSPHERES 


ELEMENTS 

IN REACTANTS, 

ATOMS 


FUEL 

OXIDIZER 

H 



2.5923 

10.34 

0 . 

C 



1.7424 

6.95 

0 . 

0 



4.7326 

0.13 

2.00 

N 



C.C526 

0.21 

C. 

ENTHALPY 

OF REACTANTS, 

KCAL 

-21.4353 




PRESSURE, ATM 

TEMPERATURE, DEG K 

HEAT CAPACITY, CAL /DEG K 

ENTHALPY, KCAL 

ENTROPY, CAL / DEG K 

FROZEN GAMMA 

SHIFTING GAMMA 

MOLS OF GAS 

MOLECULAR WEIGHT 


CHAMBER 

15.0CCC 

3441.3549 

195.8732 

-21.4353 

254.0893 

1.2128 

1.1222 

3.7754 

26.4874 


CENSITY, GM / CC C. 001406923 


100% PBAN ( 0/F = 3.11) 
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HYBRID SYSTEM NUMBER IF PRESSURE * 5 AT MB SPHERES 


ELEMENTS IN REACTANTS* 

ATOMS 

FUEL 

OXIDIZER 

H 

3.1986 

10.94 

0. 

c 

2.1499 

6.95 

0 . 

D 

4.3777 

0.13 

2.00 

N 

ENTHALPY OF REACTANTS* 

PRESSURE* ATM 
TEMPERATURE* DEG K 

0.0650 

KCAL -26.4484 

CHAMBER 

5.000C 

3264,0998 

0.21 

G. 

HEAT CAPACITY* CAL /DEG K 196.0999 
ENTHALPY* KCAL -26.4484 
ENTROPY* CAL / DEG K 277.4010 
FROZEN GAMMA 1.2244 
SHIFTING GAMMA 1.1297 
MOLS OF GAS 4.1425 
MOLECULAR WEIGHT 24.0238 
DENSITY, GM / CC C.C0C4481I2 




HYBRID 

SYSTEM NUMBER 

IF 

PRESSURE *15 

ATMS SPHERES 


ELEMENTS 

IN 

REACTANTS, 

ATOMS 


FUEL 

OXIDIZER 

H 




3.1986 

10.34 

C. 

C 




2.1499 

6.95 

0. 

0 




4.1777 

0.13 

2.00 

N 




0.0650 

0.21 

0. 

ENTHALPY 

OF 

REACTANTS* 

KCAL 

-26.4484 




PRESSURE, ATM 


CHAMBER 

15.0000 

TEMPERATURE* DEG K 


3410.4482 

HEAT CAPACITY* CAL 

/DEG 

K 166.9462 

ENTHALPY, KCAL 


-26.4484 

ENTROPY, CAL / DEG 

K 

268.3219 

FROZEN GAMMA 


1.2200 

SHIFTING GAMMA 


1.1299 

MOLS OF GAS 


4.1087 

MOLECULAR HEIGHT 


24.3386 

DENSITY* GM / CC 


C. 001304485 


100% PBAN (0/F = 2.34) 



HYBRID SYSTEM NUMBER 2A 


PRESSURE « 5 ATMOSPHERES 


ELEMENTS 

IN 

REACTANTS, 

ATOMS 


FUEL 

OXIDIZER 

BINDERS 

H 




2.2811 

9. 87 

0 . 

0 . , 

C 




1.2845 

5.96 

0 . 

0 . 

0 




4.2815 

1.50 

2.00 

0 . 

N 




0.0285 

0.22 

0 . 

0 . 

AL 




0.9448 

0 . 

0 . 

1.00 

> 

ENTHALPY 

OF 

REACTANTS, 

KCAL 

■*29.2253 




PRESSURE, ATM 
TEMPERATURE, DEG K 
HEAT CAPACITY, CAL 
ENTHALPY, KCAL 
ENTROPY, CAL / DEG 
FROZEN GAMMA 
SHIFTING GAMMA 
MOLS OF GAS 
MOLECULAR WEIGHT 
DENSITY, GM / CC 


CHAMBER 
5.00CC 
3632.5299 
/DEG K 241.7865 
>19.2253 
K 241.1667 

1.1988 
1.1277 
3.3291 
27.84C1 
C.C0C5C4732 


HYBRID SYSTEM NUMBER 2A PRESSURE =15 ATMOSPHERES 


ELEMENTS IN REACTANTS, ATOMS 


H 2.2811 
C 1.2845 
0 4.2815 
N C.C265 
AL 0.5448 
ENTHALPY OF REACTANTS, KCAL -19.2253 


FUEL 

OXIDIZER 

BINDERS 

9.87 

0 . 

0 . 

5.36 

0 . 

0 . 

1. 50 

2.00 

0 . 

0.22 

0 . 

0 . 

0 . 

0 . 

1.00 


CHAMBER 


PRESSURE, ATM 

TEMPERATURE, DEG K 

HEAT CAPACITY, CAL /DEG K 

ENTHALPY, KCAL 

ENTROPY, CAL / DEG K 

FROZEN GAMMA 

SHIFTING GAMMA 

MOLS OF GAS 

MOLECULAR WEIGHT 

DENSITY, GM / CC 


15.0CCC 
3831.5544 
207.4829 
-19.2253 
233.9875 
1.1836 
1.1221 
3.2529 
28.3947 
C. 001466416 


60% PU/40% A1 (0/F = 1.72) 
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HYBRID SYSTEM NUMBER 

2B 

PRESSURE - 5 

ATMOSPHERES 



ELEMENTS IN REACTANTS. 

ATOMS 

FUEL 

OXIDIZER 

BINDERS 

H 



2.5867 

9.97 

0 * 

0 . 

C 



1.4042 

5.16 

0 . 

0 . 

0 



9.80S3 

1.60 

2.00 

0 . 

N 



0.0314 

0.22 

0. 

0 . 

AL 



0.4523 

0 . 

0 . 

1.00 

ENTHALPY OF REACTANTS* 

KCAL 

*21.7920 







CHAMBER 




PRESSURE. ATM 



5. 00 GO 




TEMPERATURE. OEG K 



3694.7493 




HEAT CAPACITY. CAL 

/OEG K 

261.1149 




ENTHALPY, KCAL 



-22.7920 




ENTROPY, CAL / DEG 

K 


248.7089 




FROZEN GAMMA 



1.1891 




SHIFTING GAMMA 



1.1178 




MOLS OF GAS 



3.4864 




MOLECULAR WEIGHT 



26.2990 




DENSITY, GM / CC 


C* 

000472993 





HYBRID 

SYSTEM NUMBER 

2B 

PRESSURE *15 

ATMOSPHERES 



ELEMENTS 

IN 

REACTANTS, 

ATOMS 


FUEL 

OXIDIZER 

BINDERS 

H 




2.5857 

9.87 

0. 

0. 

C 




1.4042 

5.36 

0. 

0. 

0 




3.9055 

1.50 

2.00 

0. 

N 




0.0314 

0.12 

C. 

0. 

AL 




0.6523 

0. 

0. 

1.00 

ENTHALPY 

OF 

REACTANTS, 

KCAL 

-22.7920 





PRESSURE. ATM 

TEMPERATURE. DEG K 

HEAT CAPACITY. CAL /DEG K 

ENTHALPY. KCAL 

ENTROPY. CAL / DEG K 

FROZEN GAMMA 

SHIFTING GAMMA 

MOLS OF GAS 

MOLECULAR WEIGHT 

DENSITY. GM / CC C 


CHAMBER 
15. COCO 
3900.9078 
221*9965 ' 
-22.7920 
241*1806 
1*1838 
1.1123 
3.4093 
26.8422 
001374237 


60% FU/40% A1 ( 0/F = 1.28) 
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HYBRID SYSTEM NUMBER 2C 


PRESSURE - 5 ATMS SPHERES 


ELEMENTS 

IN 

REACTANTS, 

ATOMS 


PUEL 

OXIDIZER 

BINDERS 

H 




2.4351 

9.37 

0 . 

0 . 

C 




1.4310 

5.36 

0. 

0 . 

0 




4.5630 

1.50 

2.00 

0. 

N 




C.0320 

0.12 

0 . 

0 . 

AL 




0.2483 

0. 

0 . 

1.00 

ENTHALPY 

OF 

REACTANTS, 

KCAL 

*"23.2269 





PRESSURE, ATM 


CHAMBER 
5. 0000 

TEMPERATURE, DEG K 


3417.7964 

HEAT CAPACITY, CAL 

/DEG 

K 232.9202 

ENTHALPY, KCAL 


-23.2269 

ENTROPY, CAL / DEG 

K 

256.2981 

FROZEN GAMMA 


1.2C28 

SHIFTING GAMMA 


1.1171 

MOLS OF GAS 


3.6483 

MOLECULAR WEIGHT 


26.5124 

DENSITY, GM / CC 


0.000488639 


HYBRID SYSTEM NUMBER 2C PRESSURE >15 ATMOSPHERES 


ELEMENTS IN REACTANTS, ATOMS 


H 2.4351 
C 1.4310 
0 4.5630 
N C.C320 
AL C. 2483 
ENTHALPY OF REACTANTS, KCAL -23.2269 


FUEL 

OXIDIZER 

BINDERS 

9.87 

C. 

0. 

5.36 

0. 

0. 

1. 50 

2.00 

C. 

0.12 

C. 

0. 

0 . 

0. 

l.CO 


PRESSURE, ATM 
TEMPERATURE, DEG K 
HEAT CAPACITY, CAL /DEG K 
ENTHALPY, KCAL 
ENTROPY, CAL / DEG K 
FROZEN GAMMA 
SHIFTING GAMMA 
MOLS OF GAS 
MOLECULAR WEIGHT 
DENSITY, GM / CC 


CHAMBER 
15.C6CC 
3588.4941 
201.6443 
-23.227C 
248.3392 
1.198C 
1.1218 
3.5872 
26.9490 
C. 001419932 


80% PU/20% A1 (0/F = 2.00) 



HYBRID SYSTEM NUMBER 2D PRESSURE * 5 ATMOSPHERES 


ELEMENTS IN REACTANTS* 

ATOMS 

FUEL 

OXIDIZER 

BINDERS 

H 

3.2779 

9.87 

0 . 

0 . 

C 

1.7258 

5.36 

0 . 

0 . 

0 

4.2205 

1.50 

2.00 

0 . 

N 

0.0386 

0.12 

0. 

C. 

AL 

0.2965 

0 . 

0 . 

1.00 

ENTHALPY OF REACTANTS* 

KCAL ^28.0115 





CHAMBER 




PRESSURE* ATM 

5.0000 




TEMPERATURE* DEG K 

3438.2041 




HEAT CAPACITY* CAL /DEG K 227.3420 




ENTHALPY, KCAL 

-28.0115 




ENTROPY, CAL ✓ DEG K 

266.2076 




FROZEN GAMMA 

1.2C57 




SHIFTING GAMMA 

1.1204 




MOLS OF GAS 

3.9007 




MOLECULAR WEIGHT 

24.7054 




DENSITY, GM / CC 

C. 000454307 




HYBRID SYSTEM NUMBER 

2D PRESSURE =15 

ATMOSPHERES 


ELEMENTS IN REACTANTS, 

ATOMS 

FUEL 

OXIDIZER 

BINDERS 

H 

3.1779 

9.87 

C. 

0. 

C 

1.7258 

5.36 

0 . 

0 . 

0 

4.22C5 

1.50 

2.00 

0. 

N 

C.C386 

0.12 

0. 

0. 

AL 

C. 2965 

0. 

0. 

l.CO 

ENTHALPY OF REACTANTS, 

KCAL -28.C115 





CHAMBER 




PRESSURE, ATM 

15.CCCC 




TEMPERATURE, DEG K 

3608. 53C2 




HEAT CAPACITY, CAL /DEG K 193. 9320 




ENTHALPY, KCAL 

-28.0115 




ENTROPY, CAL / DEG K 

259.7420 




FROZEN GAMMA 

l • 2 C 1 0 




SHIFTING GAMMA 

1.1255 




MOLS OF GAS 

3.8357 




MOLECULAR WEIGHT 

25.1087 




DENSITY, GM / CC 

C. 001320656 





80% PU/20% A1 ( 0/F = 1.50) 



HYBRID SYSTEM NUMBER 2 E 


PRESSURE 


5 ATMOSPHERES 


ELEMENTS 

IN 

REACTANTS, 

ATOMS 


FUEL 

OXIDIZER 

H 




3.0101 

9.87 

0. 

C 




1.4347 

5.36 

0 . 

0 




4.8012 

1.50 

2.00 

N 




0.0366 

0.12 

0 . 

ENTHALPY 

OF 

REACTANTS, 

KCAL 

-26.5326 




PRESSURE, ATM 
TEMPERATURE, DEG K 
HEAT CAPACITY, CAL 
ENTHALPY, KCAL 
ENTROPY, CAL / DEG 
FROZEN GAMMA 
SHIFTING GAMMA 
MOLS OF GAS 
MOLECULAR WEIGHT 
DENSITY, GM / CC 


CHAMBER 
5.0CCC 
3256.5509 
/DEG K 226.1879 
-26.5326 
K 268.1019 

1.214C 
1. 1164 
3. 9190 
25.5169 
C.C0C477426 


HYBRID SYSTEM NUMBER 2E PRESSURE =15 ATMESPHERES 


ELEMENTS 

IN 

REACTANTS, 

ATOMS 


FUEL 

OXIDIZER 

H 




3.0101 

9.87 

0. 

C 




1.4347 

5.36 

c. 

0 




4.8012 

1.50 

2.00 

N 




0.0366 

0.22 

c. 

ENTHALPY 

OF 

REACTANTS, 

KCAL 

-26.5326 




PRESSURE, ATM 
TEMPERATURE, DEG K 
HEAT CAPACITY, CAL 
ENTHALPY, KCAL 
ENTROPY, CAL / DEG 
FROZEN GAMMA 
SHIFTING GAMMA 
MOLS OF GAS 
MOLECULAR WEIGHT 
DENSITY, GM / CC 


CHAMBER 
15.0000 
3406.9760 
/DEG K 197.2591 
-26.5327 
K 259.6C49 

1. 2094 
1.1213 
3.8645 
25.8763 
0.001388326 


100% PU (0/F = 2.27) 


112 



HYBRID SYSTEM NUMBER 2F 


PRESSURE * 5 ATMCSPHERES 


ELEMENTS IN REACTANTS, ATOMS 
H 
C 
0 
N 

ENTHALPY OF REACTANTS, KCAL 


PRESSURE, ATM 
TEMPERATURE, DEG K 
HEAT CAPACITY, CAL 
ENTHALPY, KCAL 
ENTROPY, CAL / DEG 
FROZEN GAMMA 
SHIFTING GAMMA 
MOLS OF GAS 
MOLECULAR HEIGHT 
DENSITY, GM / CC 



FUEL 

OXIDIZER 

3.6220 

9.87 

0 . 

1.9669 

5.36 

0 . 

4.5067 

1.50 

2.00 

0.0440 

31.9262 

0.22 

0. 


CHAMBER 


5.0000 
3245.2159 
/DEG K 209.7912 
-31.9242 
K 282.4401 

1.2118 
1.1209 
4.2231 
23.6792 
C. 000444591 


HYBRID SYSTEM NUMBER 2F PRESSURE «15 ATMESPHERES 


ELEMENTS IN REACTANTS, ATOMS 

H 

C 

0 

N 

ENTHALPY OF REACTANTS, KCAL 



FUEL 

OXIDIZER 

3. 6220 

9.B7 

0. 

1.9669 

5.36 

0. 

4.5067 

1.50 

2.00 

C.C440 

-31.9262 

0.12 

C. 


CHAMBER 


PRESSURE, ATM 

TEMPERATURE, DEG K 

HEAT CAPACITY, CAL /DEG K 

ENTHALPY, KCAL 

ENTROPY, CAL / DEG K 

FROZEN GAMMA 

SHIFTING GAMMA 

MOLS OF GAS 

MOLECULAR HEIGHT 

DENSITY, GM / CC C. 


15.0CCC 

3389.2556 

179.2704 

-31.9262 

273.28C4 

1.2145 

1.1266 

4.1*77 

23.9942 

CC1294C76 


100% PU (0/F = 1.73) 
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